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ABSTRACT 


VI 

VI 


In  this  report,  the  general  techniques  for  determining  the  distributions  of 
products  and  quotients  of  random  variables  are  discussed.  Some  exact  and 
asymptotic  results  pertaining  to  the  distributions  of  the  products  and  quotients 
of  certain  random  variables  which  generally  occur  as  measurement  error  are 
also  presented  together  with  their  applications.  An  extensive  bibliography  is 
included  at  the  end  of  the  report. 
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I.  GENERAL  INTRODUCTION 


In  the  applied  sciences,  problems  are  frequently  encountered  con¬ 
cerning  reliability  analyses,  measures  of  efficiency,  biometry  indices, 
etc.,  which  properly  belong  in  the  framework  of  determining  the  distri¬ 
butions  of  various  algebraic  combinations  of  random  variables.  Many 
of  the  problems  associated  with  this  general  area,  and  in  particular 
with  product  and  quotient  forms  of  these  variables,  have  been  extensively 
investigated.  Craig  ^21 j  and  others,  however,  point  out  that  these 
investigations  have  resulted. in  surprisingly  little  useful  information 
which  is  generally  available  to  the  engineer,  research  scientist,  etc. 

The  scope  of  such  problems  is  surprisingly  broad  and  investigations  into 
methods  of  analysis  date  as  far  back  as  to  K.  Pearson  [97j  ,  1910.  As  a 
collective  entity,  the  literature  concerning  product  and  quotient  forms 
of  random  variables  is  widely  scattered  and,  unfortunately,  is  devoted 
almost  exclusively  to  investigations  utilizing  very  specialized  quad¬ 
rature  methods.  As  a  result,  no  publication  has  appeared  which  presents 
a  general  theory,  methods  of  application,  and  useful  tabular  results  per¬ 
taining  to  a  variety  of  product,  and  quotient  forms  of  random  variables. 

Scattered  and  incomplete  tables  concerning  exact  distribution 
percentage  points,  approximating  distributions,  and  other  applied  topics 
are  available.  To  the  engineer  unschooled  in  the  rudiments  of  random 
variable  techniques,  these  are  often  meaningless  and  force  him  to  use 


^ In  the  body  of  this  report,  numbers  in  brackets,  £  J  ,  refer  to  the 
bibliography. 


less  desirable  deterministic  methods  or  even  perhaps  to  make  unwarranted 
assumptions  so  that  the  tabular  results  might  be  used. 

The  aim  of  this  monograph  is  to  present  a  useful  expository  manu¬ 
script  on  certain  frequently  occurring  product  and  quotient  forms  of 
random  variables  which  will  present  a  general  theory,  bring  together 
and  discuss  known  results,  and  reference  pertinent  tabular  results  so  as 
to  provide  a  useful  tool  to  the  engineer  or  research  scientist. 

In  this  respect,  the  monograph  is  devoted  to  the  following  topics: 

1)  Examples  of  applied  problems  involving  products  and  quotients.  To 
demonstrate  a  wide  range  of  applications  to  engineering  problems,  a 
number  of  varied  examples  are  discussed,  mainly  from  the  viewpoint  of 
recognizing  and  properly  posing  the  problems  in  this  general  area. 

2)  General  theoretical  models  for  determining  product  and  quotient 
distributions.  The  general  techniques  for  determining  the  distributions 
of  products  and  quotients  of  random  variables  are  presented.  This  pre¬ 
sentation  treats  the  independent  and  dependent  cases  of  continuous 
variables. 

3)  Exact  results  pertaining  to  products  and  quotients  of  random  variables 
which  generally  occur  as  measurement  error.  Considerable  attention  is 
devoted  to  the  quotient  distributions 5  normal/normal,  rectangular/normal, 
triangular/normal;  to  the  product  distribution,  normal  x  normal,  rectan¬ 
gular  x  normal,  triangular  x  normal;  and  to  certain  discrete  distri- 

1 

buttons  of  the  form  —  . 


4)  Limiting  distributions,  approximations  and  asymptotic  results.  The 

n  /  e 


general  variable  Z  = 


n 

i=l 


n  Y.  is  discussed  in  respect  to 
j-n+l  3 


fitting  the  lognormal  approximation  to  Z,  along  with  other  topics. 


5)  characterizing  properties  of  statistical  distributions.  Problems 
of  the  type:  "If  X^/X^  =  Y  and  Y  follows  a  general  F  distribution,  are 

and  necessarily  chi-square  variables?"  are  treated. 

6)  An  annotated  bibliography.  Products  and  quotients  of  random 
variables  assume  many  different  representations.  Not  only  is  one 
interested  in  ratios  and  products  of  random  variables  described  by 
specific  frequency  functions,  but  also  in  functional  forms  of  variances, 
ranges,  proportions!,  etc. 

The  bibliography  presented  herein  is  very  comprehensive  in  respect 
tc  articles  pertaining  to  the  theory  of  the  distribution  of  random 
variables  described  by  probability  density  functions.  The  remaining 
entries  in  the  bibliography  represent  a  small  representative  sample 
of  the  vast  number  of  published  articles  pertaining  to  other  forms  of 
products  and  quotients. 

The  notation  used  herein  conforms  for  the  most  part  to  what  is 
believed  to  be  a  standard  format  in  statistical  literature.  The  symbols 
denoting  random  variables  will  be  X,  Y,  Z  and  the  values  of  the  variates, 


X,  y,  z. 


Frequency  or  probability  density  functions,  p.d.f.'g  will 

be  donctod  by  small  letters  f(x),  E(y) .  and  the  respective 

cumulative  distribution  functions,  c.d.f.'s,  by  F(x) ,  G(y),  .  .  .  ,  • 

A  random  variable  obeying  a  certain  probability  density  function,  say 

p 

the  Gaussian  or  normal  probability  law,  will  be  denoted  as  N(  /*,  o r~). 
Here  }X  and  or2  respectively  denote  the  expected  value  and  the  variance 

of  the  normal  random  variable. 
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11 •  ENGINEERING  applications  involving 

QUOTIENTS  OR  PRODUCTS  OF  RANDOM  VARIABLES 


The  examples  discussed  in  this  chapter  serve  to  emphasize  the 
broad  range  of  product  and  quotient  forme  arising  in  engineering 
applications  of  statistical  distribution  theory. 

2 

2.1  Cyclic  Firing  Rate  of  the  T-l60  (20  mm)  Cannon  in  Korean  Combat 

The  need  for  establishing  the  reliability  and  functional  suitability 
of  weapons  under  combat  conditions  is  apparent.  The  following  account 
describes  a  preliminary  model  for  a  statistical  analysis  which  was  used 
to  determine  the  combat  suitability  of  the  then  newly  manufactured 
T-160  gun  installed  on  jet  fighter  aircraft  during  a  period  in  the 
Korean  conflict. 

It  was  desired  to  establish  the  cyclic  gun  firing  rate  the  T-160 
achieved  during  combat  by  assessing  sight  reticle  and  scope  camera 
film  exposed  during  missions  and  from  questionnaires  filled  out  by 
pilots  and  armament  technicians  detailing  the  number  of  rounds  expended 
on  each  sortie  from  belt  counts.  The  data  were  used  to  obtain  a  prob¬ 
abilistic  estimate  of  the  typical  cyclic  rate  of  the  T-l60  for  the 
Korean  combat  test. 

The  parameters  r  and  9  defined  below  refer  to  sortie  averages.  If 

r  =  average  cyclic  gun  rate  per  sortie  in  rounds/sec/gun, 

0  =  average  camera  speed  per  sortie  in  frames/sec„ 


George  W.  Morgenthaler ,  "Cyclic  Rate  of  the  T-l60  in  Korean  Combut," 
Combat  Film  Analysis  Newsletters,  Institute  of  Air  Weapons  Research, 
University  of  Chicago,  195^. 


.*  J  <7  H"  ■<  \  ■  *w  •’  *.  *’  >  • 


then 


(No.  of  guns  firing)  — {’—mea  r  =  (No.  of  rounds  expended).  (2-1) 


There  are,  of  course,  certain  errors  in  the  input  data  for  formula 
(2-1).  The  fi  ,'ing  frame  count  may  be  in  error  because  of  the  difficulty 
inherent  in  distinguishing  a  firing  frame  from  o  nonfiring  frame, 
particularly  with  poor  quality  film,,  The  nominal  camera  speed 
(9  =  32  fr/sec.)  is  subject  to  a  factory  setting  error,  in  addition  to 
errors  induced  by  mechanical  wear,  maintenance  irregularities,  and  to 
errors  arising  from  variations  in  temperature  and  the  camera  circuit 
voltage. 

Test  firings  indicated  that  the  average  cyclic  rate,  r,  on  a  sortie 
is  a  function  of  total  time  of  fire,  the  number  of  bursts,  the  quality 
of  the  gun,  etc.  Beihg  a  function  of  random  variables,  r  is  a  random 
variable  and  thus  would  generate  a  distribution  over  a  large  number  of 
sorties.  The  expected  value,  E(r),  is  representative  of  the  sortie 
cyclic  rates  experienced  in  Korea. 

In  the  treatment  of  this  problem,  r  was  treated  as  a  product  of  the 
dependent  random  variables  0  and  r/0  and  E(r)  was  defined  by  the  relation¬ 
ship: 


E(r)  =  E 


=  E(e)  e 


(2-2) 


where: 


I* 

coefficient  of  correlation  between  9  and  ^  , 
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V(Q)  =  variance  of  ©,  and 


v(f) 


variance  of  g 


In  this  case,  the  moments  of  the  distribution  of  r  =  (9) 

r 

were  first  obtained  by  making  the  assumption  that  0  and  g  were 
normally  distributed.  This  assumption  was  supported  by  experimental 
data.  It  was  found  that  1164  <  E(r)  <  1266  rds./min.  with  a  prob¬ 
ability  of  .90.  This  estimate  was  surprisingly  low  in  comparison 
with  the  T-160  designed  cyclic  rate  of  1500  rds/min.  It  is  apparent  that 
with  the  percentage  points  of  the  distribution  of  r  =  (0)  ^g^  ,  one 

could  examine  the  prob-  bility  of  an  E(r)  such  as  was  attained.  This 
very  low  probability  would  cause  one  to  search  for  a  plausible  ex¬ 
planation. 

In  this  case,  it  was  discovered  that  since  the  T-160  is  a  gas- 
operating,  automatic  weapon  consisting  of  barrel  and  a  rotating  drum 
with  five  chambers,  then  the  inertia  of  the  drum  would  play  an  important 
role.  At  the  outset,  the  inertia  of  the  drum  is  overcome  rather  slowly  and, 
hence,  the  gun  does  not  have  a  cyclic  rate  but  rather  an  average  cyclic 
rate  depending  upon  the  duration  of  the  burst  of  fire. 

Computation  ot  the  expected  T-160  cyclic  rate  for  Korean  type 
combat  based  on  the  Korean  T-l60  burst  length  later  proved  that  values 
of  E(r)  ="  1200  rds./min.  were  to  be  expected. 


2  .2  Selection  of  a  Space  "Workhorse"  Booster. 

Many  experts  agree  that  all  long  range  space  explorations  to  be 


attempted  in  the  intermediate  future  will  best  be  initiated  from  an 
earth  or  parking  orbit.  In  this  respect,  the  success  of  such  space 
missions  will  then  greatly  depend  upon  the  successful  completion  of  a 
very  important  logistics  operation.  This  operation  will  involve  the 
transportation  of  all  necessary  equipment  and  supplies  into  the  earth 
orbit  and  the  assembly  of  it  there. 


The  selection  of  a  suitable  "workhorse"  booster  for  this  phase  of 
the  project  represents  only  one  of  the  myriad  of  complex  decisions 
facing  the  project  planners.  But,  since  this  decision  must  be  made 
years  ahead  of  most  others  and  will  invol/e  great  expenditures  of 
money,  manpower,  and  time,  it  perhaps  represents  the  most  important 
present  day  space  industry  problem. 

Let  us  assume  that,  at  the  present  time,  the  choice  of  a  workhorse 
booster  is  restricted  to  typical  systems  such  as  the  Nova,  Saturn  C-5» 
and  the  Titan  III  systems.  The  concepts  of  these  boosters  differ  con¬ 
siderably;  and  any  operational  boosters  forthcoming  from  the  e  projects 
will  most  assuredly  differ  in  respect  to  payload  capacity,  costs, 
velocities,  reliability,  date  of  assembly  of  the  first  operational  booster, 
etc.  Thus  the  problem  of  an  "optimal"  selection  of  one  booster  system  for 
the  transportation  job  must  be  based  on  a  very  realistic  evaluation  model. 

One  meaningful  index  for  an  evaluation  model  of  this  type  is  the 

"dollar"  cost  per  pound  of  equipment  placed  in  orbit.  This  index  would 

be  obtained,  of  course,  by  dividing  the  total  system  cost,  0  ,  by 

system 

the  payload  weight  placed  in  orbit,  *pay^oa^»  or 


I  = 


w 


^system 


(2-3) 


payload 


A  closer  investigation  of  C  ,  and  W  .  reveals  that  these 

system  payload 

components  are  functions  of  variable  components.  For  example,  total 
system  costs  include  both  research  and  development  costs  and  operational 
costs  -  estimates  of  which  are  highly  uncertain.  Uncertainties  in  R  &  D 
cost  estimates  intervene  in  the  form  of  unanticipated  differences  in 
estimated  costs  of  component  parts  and  in  increased  expenditures  caused 
by  modifications  in  the  design  to  meet  new  or  revised  performance 
specifications.  Operational  costs  are  highly  susceptible  to  the  failure 
rates  of  the  component  parts  and,  more  generally,  to  the  level  of 
"sophistication"  in  the  supporting  logistics  system. 

The  weight  and  space  of  payload  placed  in  orbit  are  functionally  de¬ 
pendent  upon  various  characteristics  of  the  missile  such  as  Specific 
Impulse,  Isp,  the  weight- to-stage  weight  ratio,  and  other  aero- 


It  becomes  apparent  fr  m  these  considerations  that  equation  (2-3) 


Without  the  knowledge  of  how  to  combine  the  functional  forms  of 
the  components  of  I,  only  a  point  estimate  of  I  will  be  possible.  On 
this  basis,  the  choice  of  a  booster  system  would  be  greatly  influenced 
by  the  relative  differences  in  the  point  estimates  of  I  for  these  systems. 
On  the  other  hand,  the  c.d.f.'s  of  I,  F(I),  for  each  booster  may  be  ob¬ 
tained  with  a  knowledge  of  random  variable  techniques.  Assume,  for 
illustrative  purposes,  that  the  c.d.f.'s  of  Figure  1  are  realiatio. 

Several  additional,  equally  good,  criteria  for  the  selection  of  a 
booster  system  may  be  considered  when  the  information  of  Figure  1  is 
available.  For  example,  a  system  might  be  selected  on  the  basis  of: 

1)  the  system  expected  to  become  operational  soonest  and  for  which  I 
is  greater  than  some  specified  level  with  probability  fe  50#»  say,  or 

2)  the  system  for  which  I  is  largest  at  a  certain  prescribed  level  of  prob¬ 
ability  ,  or  3)  the  system  for  which  the  maximum  cost  estimate  is  a 
minimum,  etc. 

More  important  ,  all  of  these  criteria  may  be  evaluated  simultaneously 
with  the  aid  of  estimates  of  the  respective  c.d.f.'s.  Such  an  analysis 
is  certain  to  provide  much  more  valuable  insight  into  the  problem  of 
selecting  an  appropriate  booster  system. 

2. .3  Replenishment  of  a  Life  Support  System. 

With  the  advent  of  long-duration  manned  space  flights,  estimating 
the  changes  that  will  occur  in  the  quantities  of  certain  substances 
in  the  desired  ecological  system  has  become  more  complex. 


Ideally,  a  mathematical  model  for  general  system  analysis  of  an 


ecological  system  for  long-duration  flights  will  provide  for: 

1)  A  formulation  of  the  control  problem  from  which  optimal 
control  functions  car  he  determined; 

2)  A  preliminary  design  for  the  ecological  system  in  respect 
to  system  stability  and  control  considerations; 

3)  A  method  for  determining  the  time  required  to  restore  the 
system  to  a  suitable  balanoe  in  the  case  of  a  mishap  and 
the  time  required  initially  to  put  the  system  into 
operation; and 

4)  The  determination  of  resupply  requirements. 

3 

One  aspect  of  a  preliminary  model  designed  to  include  the  above 
considerations  has  resulted  in  the  requirement  for  evaluating  a  product 
of  random  variables.  It  may  be  stated  in  the  following  manner: 

Consider  the  amount  of  oxygen  in  a  cabin  atmosphere.  The  amount 
is  affected  by  leakage,  crew  consumption,  and  resupply  from  a  storage 
capacity.  Let: 

X(t)  =  amount  of  oxygen  (in  moles)  in  the  cabin  atmosphere  at 
time  t, 

L  =  proportion  of  rate  of  loss  of  oxygen  due  to  leakage  from 
the  cabin  atmosphere  per  time  period, 


_ 

Davis,  Henry  and  Novosad,  R.  S. ,  "Control,  Replenishment,  and 
Stability  of  Life  Support  Systems,"  Journal  of  Spacecraft  and 
Rockets,  Amer.  Inst,  of  A.  and  A.,  January  1964. 


Y1  =  rats  of  increase  in  oxygon  content  in  the  cabin 
atmosphere  from  storage  per  time  period,  and 
K  -  rate  of  decrease  in  the  oxygen  content  in  the  cabin 
atmosphere  due  to  crew  consumption  per  time  period. 


The  estimated  amount  of  oxygen  at  time  t  ia 


X(t)  -  exp  (-Lt) 


exp  (Lt) 


(Yx  -  K)  dt  +  X(0)  , 


(2-5) 


where 

X(0)  denotes  the  initial  condition. 

In  this  study,  I  -  exp  (Lt)  and  Z  =  (Y^  -  K.)  represent  random 
variables  which  are  functions  of  time.  A  knowledge  of  combining  random 
variables  in  product  forms  is  required  for  solving  this  problem. 

2.4  Stochastic  Differential  Equations  with  Product  Coefficients. 

In  the  space  industry,  stochastic  differential  equations  whose  co¬ 
efficients  are  functional  forms  of  products  and  quotients  of  random 
variables  frequently  play  important  roles.  They  appear  in  investigations 
of  control  theory  problems,  in  analyses  of  stresses  and  material  per¬ 
formances,  and  in  estimations  of  the  effects  of  hypervelocity  impact. 

Most  frequently,  the  designer  or  engineer  treats  this  type  of 
equation  as  n  deterministic  model  by  using  moments  of  the  distribution  as 
the  constant  coefficients  in  order  to  obtain  a  solution.  If  bounds  are 
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required  for  a  particular  solution,  digital  computer  simulation  of  the 
effects  of  small  perturbations  in  the  coefficients  usually  affords  a 
convenient  method  of  analysis. 

In  the  field  of  missile  trajectory  analysis,  both  very  complicated 
models,  auch  as  a  complete  design  optimization  analysis,  and  very 
simple  models,  such  as  the  path  of  motion  in  a  plane  of  a  point  particle, 
are  needed.  The  latter  prc  ?  analytic  solutions  in  closed  form  from 
which  valuable  insight  regaru  j  trajectories  can  be  obtained.  They 
also  furnish  simple  trajectory  patterns  that  are  of  great  value  for 
problems  involving  the  simultaneous  optimization  of  design  and  tra¬ 
jectory.  The  more  complicated  models  make  it  possible  to  take  into 
account  many  effects  that  must  of  necessity  be  omitted  from  the  simpler 
models  and  also  provide  a  good  framework  for  analyzing  the  complete 
system  of  simple  models. 

One  of  the  simpler  design  models  is  that  of  electrical  circuitry  analysis. 
Such  situations  arise  frequently  in  electrical  circuits  where  the  current 
flowing  in  one  circuit  is  influenced  by  the  current  flowing  in  another 
through  direct  interconnection  or  through  a  mutual  inductance.  The 
differential  equations  for  the  currents  are  obtained  by  considering  the 
potential  drops  across  vurious  elements  of  the  circuit.  This  leads  to 
the  following  pair  of  simultaneous  equations: 


h  Ti  +  Li 

I„  +  L, 


dl  M  dl2 
„  + 


=  E  cos  W  t, 
0 


dIl 

dt  +  M  -Jjr  =  0, 
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whore  the  coefficients  refer  to  standard  circuitry  notation.  In  the 
above  form  we  have  differential-integral  equations.  By  substituting 
the  definition  of  the  current,  I  =  ,  and  rearranging  terms,  we  obtain 


The  usual  method  of  analysis  of  this  model  consists  of  using  the 
expected  values  of  the  parameters  to  determine  a  solution.  The  ex- 
peoted  values  are  obtained  by  measuring  the  parameters,  such  as 
resistance,  etc.  These  parameter  values  permit  a  unique  solution  to 
the  model. 

As  the  parameter  values  obtained  in  this  manner  are  subject  to  ex¬ 
perimental  or  measurement  error,  a  far  more  realistic  model  is  obtained 
by  considering  the  parameters  as  random  variables.  A  "family  of 
solutions"  is  generated  by  considering  the  density  functions  of  the 
coefficients  of  the  above  set  of  equations.  This  family  of  solutions 
will  allow  many  additional  questions  to  bn  answered  concerning  the  prob¬ 
ability  density  function  of  the  "solution". 

2.5  Measurement  of  Badlation  by  Electronic  Counters. 


Proportional,  Geiger,  and  scintillation  counters  are  often  used  to 
detect  X  and  tS  radiation,  as  well  as  other  charged  particles  such  as 


electrons  and  oC  particles.  Design  of  these  counters  and  their 
associated  circuits  depends  to  some  extent  on  what  is  to  be  detected. 

A  device  common  to  all  counters  is  a  scaler.  This  electronic  device 
counts  pulses  produced  by  the  counter.  Once  the  number  of  pulses  over 
a  measured  period  of  time  is  known,  the  average  counting  rate  is  ob¬ 
tained  by  simple  division.  If  the  rate  of  pulses  is  too  high  for  a 
mechanical  device,  it  is  necessary  to  scale  down  the  pulses  by  a  known 

factor  before  feeding  them  to  a  mechanical  counter.  There  are  two  kinds 

4 

of  scalers  :  the  binary  scaler  in  which  the  Gcaler  factor  is  some  power 
of  2,  and  the  decade  scaler  in  which  the  scaling  factor  is  some  power 
of  10. 

A  typical  binary  scaler  lias  several  scaling  factors  ranging  from 
0  14 

2  to  2'  .  The  scaling  circuit  is  made  up  of  a  number  of  identical 

"stages'1  connected  in  aeries,  the  number  of  stages  being  equal  to  n, 
where  2n  is  the  desired  scaling  factor.  Each  stage  is  composed  of  a 
number  of  vacuum  tubes,  capacitors,  and  resistors,  connected  so  that 
only  one  pulse  of  current  ia  transmitted  for  every  two  pulses  received. 
Since  the  output  of  one  stage  is  connected  to  the  input  of  another,  this 
division  by  two  is  repeated  as  many  times  as  there  are  stages.  The  out¬ 
put  of  the  last  stage  is  connected  to  a  mechanical  counter  that  will 
register  one  count  for  every  pulse  transmitted  to  it  by  the  last  stage. 
Thus, if  N  pulses  from  a  counter  are  passed  through  a  circuit  of  n  stages  , 

V) - 

Cullity,  B.  D.  ,  Elements  of  X-Ray  Diffraction,  Addiaon-Vlesley 
Publishing  Company,  Reading,  Massachusetts,  1956. 
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only  N/2n  will  register  on  the  mechanical  counter. 

Because  arrival  of  X-ray  quanta  in  the  counter  is  random  in  time, 
the  accuracy  of  a  counting  rate  measurement  is  governed  by  the  laws  of 
probability.  Two  counts  of  the  same  X-ray  beam  for  identical  periods 
of  time  will  not  be  precisely  the  same  because  of  the  random  spacing 
between  pulses,  even  though  the  counter  and  scaler  are  functioning  per¬ 
fectly.  Clearly,  the  accuracy  of  a  rate  measurement  of  this  kind  improves 
as  the  time  of  counting  is  prolonged.  It  is  therefore-  important  to  know 
how  long  to  count  in  order  to  attain  a  specified  degree  of  accuracy. 

This  problem  is  complicated  when  additional  background  causes  contami¬ 
nation  in  the  counting  process.  This  unavoidable  background  is  due  to 
cosmic  rays  and  may  be  augmented,  particularly  in  some  laboratories,  by 
nearby  radioactive  materials. 

Suppose  we  want  to  estimate  the  diffraction  background  in  the  presence 

of  a  fairly  large  unavoidable  background.  Let  N  be  the  number  of  pulses 

counted  in  a  given  time  from  a  radiation  source;  let  be  the  number 

counted  in  the  same  time  with  the  radiation  source  removed.  The  N, 

o 

counts  are  due  to  unavoidable  background  and  (N-N^)  to  the  diffractable 
background  being  measured.  The  relative  probable  error  in  (N-N  )  is 

67 

^N-N.  =  “N^r—  P^cent. 

b  b 

Since  N  and  are  random  variables,  the  desirability  of  obtaining 
the  density  function  of  the  above  quotient  form  of  a  random  variable  is 
apparent. 


* 


2.6  The  Pert  Model. 

Project  schedules  of  many  kinds  may  be  schematically  presented  as 
a  network  of  nodes  and  connecting  arcs.  For  illustrative  purposes,  con¬ 
sider  the  simple  project  schedule  of  Figure  2. 


PERT  Representation  of  a  Simple  Project  Network. 
Figure  2. 


Each  arc,  a,  a’,  b,  .  ,  .  ,  f,  in  Figure  2  represents  an  "activity" 
which  must  be  accomplished  during  the  project  schedule.  The  length 
of  each  arc  denotes  the  time  required  to  complete  each  activity.  The 
nodes,  A,  A' ,  B,  .  .  .  ,  F  ,  represent  events  marking  the  beginning  or 
the  end  of  an  activity.  There  are,  of  course,  two  special  events  or 
nodes;  the  initial  event,  S,  and  the  terminal  event,  F,  between  which 
all  project  activity  is  contained. 

All  activities  which  must  be  accomplished  during  one  phase  of  the 
project  and  which  may  be  worked  on  simultaneously  are  shown  in  parallel. 


In  Figure  2,  activities  a  and  a'  are  such  events.  In  addition,  the 
chain  of  activities  a,  b,  and  e  may  be  worked  on  concurrently  with  the 
series  of  activities  a',  b',  d,  and  e'.  Activities  which  may  not  begin 
before  the  termination  of  others,  as  f  in  Figure  2,  are  shown  in  series. 

For  fixed  activity  durations,  a  very  simple  algorithm  gives  the 
length  of  time  required  for  the  total  project.  "PERT"  generalizes  the 
method  of  approach  to  recognize  uncertainties  in  activity  durations 
by  considering  them  as  random  variables.  The  usual  assumption  is 
that  these  durations  are  described  by  independent  distributions,  eaoh 
with  a  finite  range.  Of  critical  interest  in  the  PERT  analysis  is  the 
distribution  of  the  random  variable  describing  the  project's  duration. 

When  the  relevant  activity  durations  are  known  with  certainty, 
finding  the  project  duration  is  a  trivial  matter  even  for  very  large 
networks.  Unfortunately,  in  many  space  industry  research  and  develop¬ 
ment  projects,  the  time  durations  for  various  activities  are  known  only 
with  a  high  degree  of  uncertainty.  For  this  reason,  the  PERT  system 
was  created  to  facilitate  network  planning. 

The  basic  data  required  for  PERT  are  the  distributions  of  the 
activity  durations.  The  data  for  these  distributions  are  obtained 
from  technicians  wl  v  have  had  some  experience  with  the  type  of  activity 
involved.  The  distribution  of  the  project's  duration  is  a  composite 
of  these  activity  duration  distributions.  For  all  activities  appearing 
in  series,  the  total  activity  duration  for  that  series  is  the  sum  of 
the  random  variables  corresponding  to  each  pro ject  activity.  All  activiti 
appearing  in  parallel  are  treated  by  determining  the  distribution  of  the 


maximum  times  of  the  activities.  As  an  example,  the  simple  network 
in  Figure  2  is  treated  in  the  following  manner:  Events  A,  R,and  E 
mark  the  end  of  project  activities  a,  b,and  e.  The  total  project  duration 
time  of  these  activities  is  a  random  variable  determined  from  the  sum  of 
the  three  random  variables,  a,  b,  and  e.  Let  U  represent  this  sum  so 
that 

U  =  a  +  b  +  e. 

Activities  a',  c,  d,and  e'  are  performed  concurrently  with  the  series 
U.  Let  V  represent  this  sum  of  random  variables  so  that 

V  =  a1  +  c  +  d  +  e' . 

It  is  easily  seen  that  the  total  project  time  required  for  com¬ 
pletion  of  this  phase  of  the  project  is  a  random  variable  determined 
by  the  distribution  of  max  [».  A  ,  the  maximum  of  U  and  V. 

As  activity  f  must  begin  after  the  completion  of  all  other  events 
in  U  and  V,  the  total  project  duration,  T,  is  of  the  form 

T  =  max(U,  V)  +  f. 

Through  this  formulation,  not  only  the  distribution  of  project 
duration  may  be  investigated,  but  certain  other  equally  important 
topics  as  well.  One  such  topic  concerns  the  establish  -n.-'nt  of  starting 
dates  for  various  series  of  activities.  The  feasibility  of  starting  the 
series  of  activities,  say  U  and  V,  on  the  same  date  may  be  investigated 
by  considering  the  random  variable  Q  =  U/V.  Thi3  idea  suggests  that 
through  examination  of  such  probabilities  as 
Pr(U/V  £k)  and  Pr  (k’<  U/V  ^  k) , 
the  need  for  rescheduling  U  or  V  may  be  determined. 


III.  GENERAL  THEORETICAL  MODELS  FOR  QUOTIENTS  AND 
PRODUCTS  OF  RANDOM  VARIABLES 


The  notion  of  probability  plays  nn  important  role  in  statistical 
theory;  yet  in  a  chronological  sense,  an  adequate  definition  of  "the 
probability  of  an  event"  has  been  subjected  to  an  ensemble  of  varied 
approaches.  Therefore,  a  brief  mention  of  the  measure-theoretic  con¬ 
cept  underlying  the  theorems  and  definitions  presented  herein  is  per¬ 
haps  warranted.  The  measure  theory  approach  to  probability,  popularized 
by  Cramer  [?.2  J  and  others,  features  the  embodiment  of  the  fundamental 
notion  of  probability  in  measure-theoretic  ideas  through  the  concept  of 
a  theory  of  sets  and  the  "measure"  of  a  set.  The  theorems  in  this 
section  are  presented  with  the  implication  that  their  rigorous  formulation 
may  be  established  by  utilizing  certain  measure- theoretic  concepts. 

3-1  Cumulative  Distribution  Functions. 

In  the  univariate  case,  the  cumulative  distribution  function,  c.d.f., 
of  a  random  variable  X  is  defined  by  the  following  postulates! 

if  X1  <  x2’  then  F(xo)  "  F(x  )  >0,  (}_1) 

F(-oo)  =  u,  F(+oo)  =  1,  in  the  limit  sense,  (3-2) 

limit  F(b)  =  F(x)  .  (3-3) 

The  notation  of  (3-2)  implies  that  the  limit  F(x)  exists  as 
x  *•  (-oo)  or  (+oo).  Since  (3-1)  defines  F(x)  as  being  monotonic  ,it 
follows  that  F(x)  has  at  most  an  enumerable  set  of  discontinuities  and 


that  the  limits  F(x+)  and  F(x-)  exist  everywhere.  The  values  of  F(x) 
at  discontinuities  are  fixed  by  (3-3) »  It  follows  from  (3-1)  and  (3-2) 
that  F(x)  is  non-negative. 

The  relation  between  the  probability  statements  about  the  random 
variable  X  and  its  c.d.f.  is  expressed  by 

Pr(X  £  x)  «  F(x).  (3-4) 

Two  important  classes  of  o.d.f.'s  may  now  be  characterized; 

i)  a  discontinuous  c.d.f.,  F(x) ,  oharaoterizes  a  relation  such 

that  each  member  x^  of  an  at  most  enumerable  set  of  points 

x^,  x^,  ...  is  associated  with  a  respective  probability 

p,  2  0,  such  that  ^  p.  =  1 ,  and  that,  the  follov:ing  condition 
i  1 

holds; 

Pj»  for  Xj x.  (3-5) 

d 

ii)  the  second  important  class  of  c.d.f.'s  is  characterized  by 
the  existence  of  the  function  f(x)  5:  0  such  that 


Pr(X<  x)  =  F(x) 
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2 .1.2  The  Bivariate  Case:  If  R  denotes  a  rectangular  area  in  the 


2 

(x.  y)  plane,  x  <  x2’  ^i  y  <  Vg  i and  F(x,  y)  denotes  the 

second  difference, 


,a(x,}')  -  F(x2,y2)  +  F(x1,y;i)  -  FU^y^  -  Fix^y^), 


(3-8) 


then  the  joint  c.d.f. ,  F(x,  y) ,  of  random  variables  X  and  Y  is  sub¬ 
jected  to  the  following  postulates! 

^  F(x,  y)  SO  (3-9) 

and 


F(-oo ,  y)  a  F(x,  -oo )  =  0,  F( +oo ,  +oo  )  =  1.  (3-10) 

By  allowing  x^,  y^  - ►  -oo  in  (3-9) »  we  may  conclude, using 

(3-10), that 

F(xn  y2)  -  F(x1$  yx)  >  0,  if  y2  =*  yx  .  (3-11) 


Similarly, 

F(x2,  y^)  -  F(x1,  y  )  >0,  if  x2=>  x^ 

From  these  postulates,  it  is  concluded  that  F(x,  y)  is  monotonio  in 
each  variable  and  the  limits  F(x  +,  yj,  F(x,  y+)  exist  everywhere.  It 
is  easily  shown  that  F(x,  y)  is  discontinuous  on  at  most  an  enumerable 
sot  of  lines,  x  =  constant,  (similar  results  for  y) .  Finally, if 

X1  - *  ~°°  *  an(*  y±  - *"  ^-n  (3-9),  then  F(x,  y)  >  0  because  of 

(3-10).  The  values  of  F(x,  y)  at  the  discontinuities  are  fixed  by  (3-12) 
so  that 
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F(x,  y)  =  F(x  +,  y)  =  F(x,  y  +) 


012) 


Again  the  connection  between  the  probability  statements  about 
random  variables  X  and  Y  and  their  joint  c,d.f.  is  determined  by 

Pr(X  <x,  Y<y)  =  F(x,  y).  (3-13) 

In  the  bivariate  oase,  again  the  discontinuous  and  continuous  cases 
are  of  particular  interest^ 


1.  The  discontinuous  case  is  characterized  by  the  existence  of 


an  at  most  enumerable  set  of  points  (x^,  y^),  1  =  1,  2,  .  .  .  , 

and  associated  probabilities  p.  such  that  ^  p ,  =1,  So  that  F(x,y)  is 

l  i  1 


F(x,y)  Pj, 


X.SX, 


yj  -  y 


(3-14) 


2.  The  continuous  case  implies  that  there  exists  f(x,y)  >  0 
such  that 


F(x,y)  =  If  f(\  i  *>2>  dv 


d  ^ 


-00  -oo 


and 


Pr(X,  Y  C  fi) 


’/  fKx' 


y)  dxdy  . 
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(3-15) 


(3-16) 


The  mixed  case  which  is  treated  through  application  of  the  Stieltjes 
integral  will  not  be  covered. 
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3  .2  Marginal  Distributions. 

The  marginal  distribution  of  X  associated  with  the  joint  c.d.f. 
of  random  variables  X  and  Y  is  defined  by  the  relation 

Pr(X«x)  =  P(X  <  x,  Y  «=  +oo )  =  F(x,  +00).  (3-1?) 

Hence  F(x,  +00)  is  the  univariate  c.d.f,  of  X  and  is  called  the 
marginal  distribution  function  of  X.  Similarly  F(+oo,  y)  is  called 
the  marginal  distribution  function  of  Y. 

For  the  discontinuous  joint  c.d.f., 


F(x  ,+00)  =  ^  ^  p.. 


X, 


(3-18) 


and  for  the  oontihuous  case 


F(xv+oo)  a  I  J  to,v  V4l»idi2-/’ 


-00  ~oo 


-00 


(3-19) 


3 '3  Statistical  Independence. 

If  F(x,  y)  is  the  joint  c.d.f.  of  X  and  Y,  then  F(x)  =  F(x,+oo) 
and  F(y)  =  F(+oo  ,  y)  as  already  noted.  The  random  variables  X  and  Y 
are  said  to  be  statistically  independent  if  and  only  if 

F(x,  y)  =  F(x)  F(y)  (3-20) 

which  in  turn  implies  that  the  following  relations  hold: 


f(x,  y)  =  f(x)  g(y)  for  the  continuous  case 


(3-21) 


P(X  a  x,  Y  =  y)  =  Pr(X  =  x)  Pr(Y  =  y) 


(3-22) 


for  the  discrete  case. 

Ab  a  result,  the  following  definition  may  be  made: 

Definition:  The  random  variab3.es  X  and  Y  are  stochastically  inde¬ 
pendent  if  and  only  if  f(x,  y)  5  ,f(x)  g(y).  If  f(x,  y)  oannot  be  ex¬ 
pressed  identically  as  the  produot  of  the  marginal  distributions,  then 
X  and  Y  are  not  statistically  independent. 

3  .4  Transformation  of  Variables. 

Change-of-variable  integration  techniques  are  used  frequently  in  the 
Btudy  of  algebraic  combinations  of  random  variables.  In  the  following 
discussion,  only  the  continuous  random  variable  case  is  treated  as  the 
discrete  case  is  analogous  and  presents  no  additional  difficulties.  At 
this  point,  it  may  be  helpful  to  present  a  special  problem  to  emphasize 
the  ideas  involved. 

Consider  a  random  variable  X  described  by  the  p.d.f., 
f(x)  =  M-x,  0  «=  x  «=  2  , 

=  0,  elsewhere. 

The  random  variable  X  is  defined  on  a  set  =  j  0  «=  x  ■=  2  }  where 

f(x)  s-O.  Define  the  random  variable  Y  as  an  algebraic  combination  of 

13  13 

X,  say  Y  =  X  ,  and  consider  the  transformation  y  =  yr  .  Under  this 

transformation  the  set  is  mapped  into  the  3et  $  -  j  0  «*=  y  gr  }  ; 
moreover,  the  transformation  is  one-to-one.  The  one-to-one  corres¬ 
pondence  between  the  points  of  and  (3  insures  that  for  every  event 
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0  «  a  «<=  b  *5  J  ,  the  event  a  c  Y  «=  b  will  occur  when  and  only  when 
the  event  kyf~n  cX  <  k\j  b"  occurs. 


Thus, 


3, -  3 _ 

Pr(a  Y  b)  -  Pr(  4\  a  *=  X  4  <n 


■/ 


43VT 


)i  xdx  . 


By  changing  the  variable  of  integration  in  the  above  so  that 


x=  4^ 


dx  _  4 

Jty  “  “  ^ following  result  ia  obtained: 


Pr(a  <?<b) 


7T73)  d]r 


— 73  <*y 

3  y° 


The  p.d.f.  of  Y  is 


g(j)  =  ^-731  0  *r  y  -e  | , 


0,  elsewhere. 


The  change  of  variable  technique  in  the  univariate  case  is  summarized 
by  the  following  theorem: 


.Theorem  1:  Let  X  be  a  continuous  random  variable  described  by  the  p.d.f., 


f(x),  defined  on  the  set  V  =  |  a  <x  <b|  and  let  y  =  Mx)  be  a 
monotonia  transformation  having  a  unique  inverse  x  =  P  (y)  so  that 
under  the  transformation,  the  set  V  maps  by  one-to-one  correspondence 
into  /9  =  [  a'  £  y  r^b'  J  •  Further,  let  (x)  exist;  then  the  p.u.f. 

of  Y  is  given  by 


g(y)  =  fx  (  ^"]'(y)) 


£-1(y)  .  |a' s  y 


(3-23) 


5.4.1  Bivariate  Case;  The  method  of  finding  the  p.d.f.  of  one  vari¬ 
able  of  a  continuous  type  may  be  extended  to  a  function  of  n  random 
variables.  Two-variable  transformations  are  considered  now  as  a 
special  case. 

Allow  y^  =  h(x^,  x^)  and  y^  =  g(x^,  x^)  to  be  a  one-to-one  trans¬ 
formation  of  random  variables  and  with  existing  continuous  first 
partial  derivatives.  This  transformation  will  map  a  two-dimensional 
set  )S  defined  in  the  (x^,  x^)  plane  into  a  two-dimensional  set  (3 
defined  in  the  (y^,  y^)  plane.  Under  these  conditions  the  following 
theorem  holdsJ 

Theorem  2:  Let  X^  and  be  continuous  random  variables  with  joint 
p.d.f.  f(x^,  Xp)  defined  on  the  two-dimensional  set  V  =  |a^x^^b,  c^x^d}  . 
Let  y.^  =  hCx^,  x^)  and  y^  =  g(x^,  Xp)  be  a  transformation  with  a  unique 
inverse  x1  =  h  1(y1>  y and  x^  =  g  1(y1,  y?)  so  that  the  first  partial 
derivatives  of  the  inverse  functions  exist  und  the  set  2>  is  mapped 
under  the  transformation  into  ft  -  {n'f  y^ib',  c'  iy0id'|  In  a 
one-to-one  correspondence.  Then  the  joint  p.d.f.  of  random  variables 


Yx  and  Y2  is 


28 


Frequently  in  the  treatment  of  the  two-dimensional  case,  the  random 
variables  1  are  treated  as  being  independent  so  that  f  (  •  ) 

becomes  fx  (h  1(y1,y2)).  fx  (g~1(y1 ,y2) ) .  In  addition,  the  transformation 

function  Y.^  is  usually  set  equal  to  the  algebraic  combination  of  in¬ 
terest,  say  a  quotient,  so  that  Y^  =  The  random  variable  Y^  is 

defined  as  Y.,  =  X^.  The  p.d.f.  of  Y^  =  X^/X2*  the  quotient  of  interest, 
is  then  obtained  by  determining  the  marginal  p.d.f.  of  Y^.  This  pro¬ 
cedure  in  most  cases  greatly  simplifies  the  calculations  involved  in 
finding  a  desired  p.d.f. 

Theorem  2  is  essentially  a  corollary  to  a  more  general  theorem 

dealing  in  several  variables.  This  general  theorem  is  summarized  by 

stating  that  if  a  set  of  transformation  functions,  y,  »  u. (x, ,. . . ,x  ) , 

i  i  1  n 

i  --  l,...,n,  exists  with  inverse  functions,  xi  a  . . .  ,y  ),  and  if 

the  conditions  of  one-to-one  transformations  and  the  existence  of  first 
partial  derivatives  hold*,  then 


where 
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Aa  an  example  of  the  use  of  (3-24),  we  shall  consider  the  following 
textbook  type  example!  2  2t 

“VAX-  4*  X_  ) 

1  Id 

Allow:  f(x1,x2)  dx^dx^  =  ™  e  dx1<lx2  and 


conaid'ir  the  transformation: 


-1  x2 

*2  "  . 


The  domain  of  the  transformation  is  y,  2:  0,  0  <  y  -c  2  If  .  The  in- 

X  u 

verse  of  the  transformation  is 


Thus , 


3 

(nos  y2)(-y1  sin  y2> 

1  d 

(sin  y2)(y1  cos  y2) 

hence,  from  0-24),  the  joint  probability  density  function  ofjY^, 
is 

1  yl2 

f(ylty2)  -  r-^r  o  y1  dy2> 


3  Transformations  Not  One-to-One;  Let  f(x^,*2,  .  .  .  ,  x^)  be  the 

n-dimensionul  p.d.f.  of  continuous  random  variables  X..  ,  X.,  .  .  .  ,  X  • 

i  d  n 

Under  the  transformation  y^  =  u^Cx.^,  Xg,  .  .  .  ,  xr) , 

y2  =  u2Cxr  x2 . Xn} . yn  un(xl’  X2 . V  ’  tho 

n-dimensional  spaoe  where  f(x^,  x^,  •  .  .  ,  x^)  >  0  is  mapped  into 
ft  in  the  (y  ,  y2,  .  .  .  ,  yQ)  space.  Under  transformations  whioh  are 

not  one-to-one,  each  point  of  V  will  correspond  to  one  point  in  ft  , 
but  to  certain  points  in  /Q  there  will  correspond  more  than  one  point 
in  if  .  The  difficulty  presented  by  this  circumstance  is  diminished  in 
the  following  manner. 

If  the  set  %  may  be  represented  as  a  union  of  r  mutually  exclusive 
sets  .  .  .  ,  so  that  the  transformation  y^,  y^,  .  .  .  ,  y^ 

defines  a  one-to-one  transformation  of  each  subset  ^  into/?  .  Then 
the  groups, 


X1  =  Uli  (yl’  y2’  *  ’  *  ♦  yn)f 

Xg  -  (y^,  y^,  ...  i  y^) ,  (i  =  1,  2,  •  .  •  ,  r) ,  (3-26) 

*n  *  "ni  1'  y2 . >r»)' 
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represent  r  groups  of  n  inverse  functions.  If  the  first  partial 
derivatives  exist  and  are  continuous ,and  if 


II 

•H 

*  uii 

3  uli 
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3y2 

9 

yn 

Bu2i 

3  U2i 

a 

U2i 

9yl 

3  y2 

9 

yn 

(i  =  1 ,2 ,  •  .  .  ,  r), 
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are  not  identically  zero,  the  ehange-of- variable  technique  may  be 
applied  to  the  union  of  the  r  mutually  exclusive  subsets  of  V.  The 
n -dimensional  p.d.f.  of  random  variables  Y^,  Y2  •  •  •  ,  Yq  is  then 
given  by 


8(y1iy2»«*«iyn)  K±l 

i»i 


^1  ’  " 


,X  <ull 


•V . Vyi’ 


The  marginal  p.d.f.  of  any  one  Y,  say  Y^,  becomes 

g(y2<  •  •  •  i  yQ)  dy^  .  .  .  dyn  . 

A  frequently  used  example  will  serve  to  illustrate  this  result.  Let  X 
follow  a  norma],  p.d.f.  with  parameters  (p  =  0,  <r“  =  1 )  .  Consider  the 


*1(»1) 


/CD  /•  0 

‘  '  7 

-UD  -00 


This  is  accomplished  by  observing  that  in  the  case  of  a  continuous 
random  variable,  the  Pr(X  =  b)  =  0;  we  may  define  the  p.d.f.  of  X  at 
any  point  b,  or  for  that  matter  any  set  of  points  o C  with  measure  zero, 
i.e.,  the  property  that  Pr(X  C  oC)  =  0,  without  affecting  the  distri¬ 
bution  of  X.  Accordingly,  two  mutually  exclusive  sets,  and 
are  obtained  by  defining  the  normal  p.d.f,  ac 
1  1/,  2 

f  (x)  =  ■  jin—  e  X  ,  -oo  «e  x  0,  0  «=  x  c  oo, 

V  27r 

=  0  ,  x  =  0  , 


The  function  y  =  x  now  defines  a  one-to-one  transformation  which 
maps  each  ^  and  into  |  0  -=y  -=  oo}.  Tne  two  inverse  functions 

of  y  =  x2  are  x  =  -  ^"y  and  x  =  ^  y  so  that  =  -X£  and 

Using  (.3-27),  we  have  the  well  known  chi-square  distribution 


with  one  degree  of  freedom, 


.1 

2 

e  ,  0  c  y  -=  oo 

The  fundamental  theorems  involving  the  algebraic  product,  quotient, 
sum,  and  the  difference  of  random  variables  X^  and  may  now  be  stated. 

These  theorems  are  a  consequence  of  (3_24)  and  (3-2?).  (See  footnote  8,  p.  42. 

Theorem  3:  The  random  variables  product  -  X^  will  be  distrib¬ 
uted  as 

^2  ■  f\,X,  (Vy2’  *2> 

e  1 

s  0, elsewhere 

provided  the  random  variables  and  are  distributed  in  accordance 
with  the  joint  p.d.f.  fY  Y  (x.  ,  x_)  which  has  been  defined  on  the 

a^ )  A2  •*. 

set  H  so  that  the  transformation  y-j  -  u^(x^,  x^)  and  -•  u^(x^,  X2) 

maps  the  set  into  the  set  in  a  one-to-one  cori'espondence. 

A  very  important  result  pertaining  to  the  product  distribution 

of  two  independent  normally  distributed  variables  may  be  obtained  by  using 

2  2 

Theorem  J>.  Assume  X  ^  is  N(0,  0^  )  and  X2  is  N(0,  cc,  )  and  defined  so 
that  PrCX^  =  0)  =0.  The  set  =|-oo  «=  00  ;  -00  «=  X2  «=  0,  0  **=  X2 -=  00 j 

is  mapped  into  the  set  p  -  |  -00  «=  y  ^  •=  00  ;  -00  ^  y  2^  0-=y2«=roo] 

in  a  one-to-one  correspondence  under  the  transformation  y^  =  *2,  y 2  =  x^m 


1_ 

*2 


dy 


2’ 


(3-28) 
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By  Theorem  3,  the  p.d.f.  of  the  product  Y^  =  X ^  is 


)  Bf - i 

J  2ir  a: 


i  *2  I  y2 


/  yl2  y2 

exp  -fcf  -2—2  +  ~~2~ 

\y2  **1  2 


Since  is  an  even  function  in  y^,  ^(y^)  “ay  be  written  as 


^(y  )  a  - 1 -  / 

7T  o:  cr„  */ 


\y2  °i 


2  «-2 
<31 


°i°2  * 


The  following  substitution  simplifies  the  above  expression. 


(# 


then  jZ’(y^)  becomes 


,>=-*-  / 
x  0  7T  re'  *r  J 


CD  2 


1 —  llyl^ 

°i  °k  \  u 


+  u  I 


!r1  «2  0 


Hence,  substituting  0=6^,  ^(y^)  may  be  written  as 


sit 


Watson,  G.  N.,  Theory  of  Bessel  Functions,  Cambridge  University 
Press,  Cambridge,  1922,  p.  181.  ’ 
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oo 


«*!>  = 


1 

2  t  <rL 


•5 


e 


«itr2 


2*°i  «2 


iyll 

-  -— —  coaht 
°i  °2  dt 


_-l_K  (— )' 

T«i<5  \*i  °2  / 


yi  ^  o, 


(3-29) 


where  K  (*)  is  a  modified  Bessel  function  of  the  second  kind  of  zero  order, 
o 

This  result  is  well  known  and  is  discussed  in  considerable  detail 
in  Chapter  IV.  A  tabulation  of  K0(*)  is  presented  in  |l4lj  and  [l42j  . 

Theorem  4;  The  random  variable  quotient  Y^  *  k^/k^  will  be  distributed  as 


«yl> 


J  \,tz  <yl  V  y2J  |y2'  *2 
0 


(3-39) 


provided  that  the  random  variables  X^  and  Xg  are  distributed  in  accord¬ 
ance  with  the  joint  p.d.f.,  fY  Y  (x1t  *  )  which  is  defined  on  the 

y  JL  c 

set  so  that  the  transformation  y^  *  u^(x^,  Xg) ,  =  Ug(x^,  Xg)  maps 

the  set  into  the  set  0  in  a  one-to-one  correspondence. 

Theorem  4  may  be  used  to  derive  the  p.d.f.  of  the  ratio  of  two 

7 

correlated  chi  variates-.  Consider  a  bivariate  normal  p.d.f, t  in 


Several  properties  of  the  bivariate  clu  p.d.f.,  have  been  extensively 
investigated  by  P.  R.  Krishnaiah,  et  al.  [69]  ,  [70]  •  N.  L.  Johnson 
[62]  and  D.  J.  Finney  [34]  present  certain  of  the  many  important 
applications  of  this  statistic. 
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random  variables  X.  and  Xp  with  parameters  (/i  =/i  =  0 ,  <T  = 

cT  =  o^,  f  ).  S.  Bose  [12]  has  shown  that  the  statistics 


and  = 


where : 


D  =  23  (xu  '  xi)2/  ®i 


2  '  V  *  £  <X2j  '  X2V  “f 


are  jointly  distributed  in  what  is  known  as  the  central  bivariate  correlated  chi 
distribution.  P.  Krishnaiah,  et  al.  [70]  ,  beginning  with  Bose's 
expression,  have  shown  that  the  joint  p.d.f.  of  U  and  V  is 


lT(n/2)+il  ^2i(U  V)  n+2l~1  exp  [ -(U2+V2)/2  (1-  ?2) 


o  ^  1  1  mvc/txi  w  VU  v/  -\u  *rv  //  c.  \x-  ^  / 

f (U,V)  =  4  (1-  ^2)n/2  y  — L - J - L - = 

i=0  r  (n/2)i !  |2(n/2)fi  r[^(n/2)+i]  (1-  ^2)n/2+i  j-  2 


By  Theorem  4,  the  marginal  p.d.f.  of  =  U/V,  under  the  transformation 
y^  =  U/V,  y2  =  V,  may  be  derived  from  (3-31). 


o©  r-» [ t  /-A.-il^i  ,  -.n+21-1.  sn+2i-l 

2  «/?r>  1  (n/2)+i  ?  (y,y?)  (y?) 

-  Ml- O  2  r(ii/2)  i,  1 2(V2)+i  r[-(„/a)ti  j  a.^7prp 


exp  K(y1y2  )2  +  y  2)/2(l  -  f*)|  jy  \  . 
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33 


Theorem  j?:  If  random  variable  X  ia  stochastically  distributed  in 


accordance  with  the  p.d.f.  f(x)  and  Y  in  accordance  with  the  p.d.f. 
g(y),  then  the  random  variable  sum  S  =  X  +  Y  will  be  distributed  in 
accordance  with 


oo 

^(s)dx 


1,  where  ^(s) 


f(s-y)  g(y)  dy. 


(3-32) 


The  well  known  fact  that  the  sum  of  two  normally  distributed  random 

variables  is  also  normally  distributed  is  shown  by  using  Theorem  Let 
2  2 

X  =  N(  u  ,  c r  )  and  Y  =  N(  p.  ,  er  ).  The  random  variable  S  -  X  +  Y  is 

xx  y  y 

distributed  in  accordance  with  (3“32)  by  Theorem  5»  Thus 


e  dx.  (3-33) 


It  follows  with  some  manipulation  that 


g(e)  = 


°2 


c 

/ 


-oc 


00 


-*( 


X  - 


3^) 
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The  expression  in  the  braces  in  (3-3*0  is  equal  to  1  by  the  properties 

of  the  normal  distribution.  Therefore,  it  follows  that  g(s)  is  normally 

p  2  2 

distributed  with  parameters  ~  p.  +  u  and  iT  +  <T  =  d"s 

b  x  y  x  y 

Theorem  6;  If  random  variable  X  is  distributed  in  accordance  with 

the  p.  d.  f.  f(x)  and  Y  in  accordance  with  the  p.  d.  f.  g(y)  and  if  X  and  Y 


are  statistically  Independent,  then  the  random  variable  difference 
W  =  X-Y  is  distributed  as 


/oo 

*<« 


)dw  =  1  , 


where 


/oo 

f(w+y)g(y)dy. 


In  Theorems  3-6,  the  fact  that  X^  and  Xg  are  assumed  to  be 
statistically  independent  allows  a  convenient  integral  expression  for 
each  algebraic  combination  of  X^  and  X^  to  be  obtained.  The  difficulties 


are  compounded  when  this  assumption  is  improper.  In  this  case,  one  is 
faced  with  the  necessity  of  obtaining  a  holomorphic  expression  for  the 
joint  density  function  of  the  two  random  variables. 

This  difficulty  may  be  illustrated  by  considering  the  following  ox- 

p  p 

ample:  Let  random  variable  X^  =  N(  ). x x  -  ;u^,  cr  =  or  )  and 

2  2 

Xp  =  N(  /x  =  cr  =  )  be  stochastically  dependent.  ThiB  de- 

pendence  is  character! zied  by  the  coefficient  of  correlation,  ^  .  The 


joint  p.d.f,  of  X1$  X 2  may  be  expressed  as 


^»(x1,x2)  = 


2  7T(r 


1 

1 

1  1  l=»g  ' 

V1-*  1 

2(l-p2) 

\  «i  r?v 

Suppose  the  p.d.f.  of  the  random  variable  quotient  Y^  =  X^/X^, 

where  Pr(X2  =  0)  =  0,  is  desired.  We  may  utilize  the  methods  of 
Section  3. 2. 4  by  defining  the  transformation, 

yl  =  X1^X2  1 
y2  15  X2* 

The  IJI  =  |y2|  and  by  Theorem  2,  the  joint  p.d.f.  of  (y^y  )  is 

ly2  I 


^(y1.y2)  = 


2T 


°i 


exp 


2(l-p)‘ 


fyi  y2  -  ^1 

‘’L 


By  integrating  (3-?5)  with  respect  to  y,,  the  marginal  p.d.f.  of  the 
random  variable  Y1  =  X.j/X2  is  obtained.  After  much  detail, the  following 
cumbersome  result  is  obtained} 


$ 


1 
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S' 

VI1 
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?<*!> 


i_  5  fi  fZ 
T  k2  rf-ztWz  +  si2) 


exp  - 


i  Ai  -  l?iv*L+  ^kl\  + 

2(l-p2)U2  ^  of) 


..  W~y i  ^2r 

+  exp{  -)?  - r - ^ — 

/ «-  *- 


^  ^2  6"l)+yl<r2(^  Vl"  V2} 


(<3i2  yi2  “2eyi  +  °iM  7r(<522  yi2  “2  e  yi  ^  °2 +  °i  > 


2,  3/ 2  * 


/ 


<rx  <g  -  ^2  V  +  yi  ^2  °i  ~  fl  (r2) 

<rx  cr2  {(i-^Kcr/  y/  -  2  P  y3.  °i  °i  + 


-Va  u  du , 


036) 


where  the  integral  function  is  the  error  Integral. 

This  example  serves  to  illustrate  that  although  in  theory  the 
change-of-variable  integration  techniques  are  always  applicable;  they 
often  result  in  very  unwieldy  integration  problems.  In  the  next  chapter, 
then,  some  of  the  more  frequently  occurring  forms  of  products  and  quo¬ 
tients  of  random  variables  will  be  investigated.  Ways  and 
methods  of  obtaining  the  desired  marginal  distributions  which  avoid 
certain  integration  difficulties  will  be  investigated. 


“Theorems- 3-b  are  widely  known  and  are  frequently  quoted 
in  statistical  literature,  no  one  person  is  credited  with 
their  derivation.  However,  Huntington  [60]  presents  a 
geometric  proof  for  each  of  these  theorems  in  his  paper. 
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IV*  PRODUCT  AND  QUOTIENT  FORMS  OF  MEASUREMENT  ERROR . 


4.1  Measurement  Error 

Physical  und  economic  limitations,  the  ever  present  constraints  in 
any  industrial  environment,  often  require  that  the  formulation  of  un 
estimate  of  a  ratio  or  of  a  product  of  two  quantities  be  based  on 
individual  measurements  of  the  two  components.  It  is  commoh,  for  ex¬ 
ample,  to  determine  the  proportion  of  a  reactant  which  has  reacted  in  a 
chemical  process  by  measuring  the  residue  after  the  process  has  com¬ 
pleted,  or  to  estimate  the  efficiency  of  a  rocket  engine  by  determining 
the  ratio  of  fuel  energy  input  to  thrust  output.  Estimates  of  efficiency, 
velocity,  and  proportions  are  examples  of  indices  which  are  necessarily 
ratios.  Product  forms  frequently  arise  in  assessing  the  probability  of 
a  successful  event,  such  as  the  launching  of  a  missile,  which  depends 
on  the  successful  operation  of  several  independent  component  subevents. 

As  is  frequently  the  case,  further  complications  arise  when  the 
individual  measurement  of  each  component  of  the  product  or  ratio  is  sub¬ 
ject  to  error.  This  error  most  often  is  due  either  to  what  in  reality 
are  true  random  fluctuations  in  the  measurable  quantity  or  to  instrument 
error  caused  by  the  inherent  limitations  of  an  imperfect  measuring  device. 
These  two  sources  of  error  are  often  treated  alike  statistically  arid  are 
colled  measurement  error.  When  the  two  quantities  are  subject  to  measure¬ 
ment  error,  their  product  or  quotient  is  distributed  about  its  expected 
value.  In  the  following,  two  examples  are  presented  which  demonstrate 
two  commonly  occurring  situations  which  require  the  use  of  approximate 
methods  to  establish  the  precision  of  an  estlmute  of  measurement  error. 


•n  tn 

m 


Whenever  economically  feasible,  the  distribution  of  error  of  each  com¬ 
ponent  quantity  may  be  determined  from  repeated  experimentation.  On  other 
occasions,  repeated  trials  are  impossible  or  undesirable?  whereupon  ,the 
precision  of  a  specific  trial  must  be  established  by  theoretical  arguments. 
This  requirement  is  frequently  necessary  in  the  acceptance  trials  of  certain 
types  of  heavy  industrial  machinery  such  as  steam  boilers,  petroleum  pro¬ 
cessing  equipment,  blast  furnaces,  generators,  large  rocket  engines, 
etc.  Such  trials  are  lengthy  and  expensive  processes  which  usually  require 
that  an  estimate  of  the  efficiency  of  the  equipment  be  made  from  a  single 
trial  operation.  It  is  necessary  from  the  point  of  establishing  confidence 
limits  to  obtain  the  probability  density  function  of  the  estimate  of 
efficiency.  Approximate  methods,  such  as  the  assumption  of  normality, must 


suffice  for  this  purpose  since  more  desirable  exact  methods  are  clearly  un¬ 
attainable. 

The  second  situation  Involves  instrument  error.  Measurements  taken 
with  several  different  measuring  devices,  all  of  a  specific  type,  generally 
are  different.  Repeated  measurement  taken  with  any  one  of  these  instru¬ 
ments  is  physically  possible,  yet  is  often  pointless;  the  instrument  has 
been  manufactured  to  certain  tolerances  and  is  incapable  of  recording  with 
a  greater  or  lesser  degree  of  precision  on  any  successive  trial.  Other 
measuring  devices  will  record  with  the  same  degree  of  precision  on  any 
specified  trial  but  in  most  cases  slight  differences  will  be  noted  in  tne 
expected  values  and  dispersion  of  measurements  taken  with  the  individual 
instruments.  The  experimenter  will  usually  have  at  his  disposal  only  one 
or  two  measurement  instruments  of  a  specified  type;  thus  he  will  be  unable 
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to  conduct  a  large  number  of  experiments  with  which  to  calibrate  his  own 
equipment.  As  a  result,  the  tolerance  associated  with  a  measurement 
accomplished  with  a  "typical"  measuring  instrument  is  best  considered  as 
being  random  in  nature  and  should  be  treated  from  a  statistical  point  of 
view. 

.In  many  industrial  applications,  the  errors  of  measurement  are  con¬ 
sidered  to  be  described  by  one  of  three  p.d.f.'s:  the  rectangular,  the 
triangular, or  the  normal  density  function. 

Instrument  error  is  often  approximated  by  the  rectangular  p.d.f. 
Certain  measurements  such  as  length,  time, and  weight  are  usually  recorded 
in  terms  of  deviations  from  some  preselected  value,  X'.  These  deviations 
y  =  (x^-  :,c')  are  assumed  to  be  in  an  interval  S  =  |a<  y  <  b|  so  that  all 
deviations  are  equally  likely.  It  is  implied  that  the  deviation  between 
the  true  value  of  the  measurement  and  the  preselected  value  is  equally 
likely  to  be  anywhere  in  the  interval  S.  A  measurement  which  has  been 
rounded  off  from  a  more  precise  measurement  is  an  example  of  an  estimate 
which  is  subject  to  rectangular  error. 

The  triangular  p.d.f.  describes  measurement  error  resulting  from 
summing  or  taking  the  difference  of  two  readings  which  are  subject  to 
rectangular  measurement  error. 

The  p.d.f.  describing  sampling  error  and  the  distribution  of  the  ex¬ 
pected  value  of  a  sample  of  measurements  is  often  taken  to  be  normal  or 
Gaussian.  Product  and  quotient  forms  of  these  three  common  p.d.f.'s  de¬ 
scribing  measurement  error  will  be  discussed  in  the  following  sections. 
These  three  r.v.'s  will  be  denoted  in  an  obvious  notation  as  R,  T,  and  N 
hereafter. 
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4.2  The  Normal  Density  Function. 


•'i.2.1  The  Product  of  Two  Normally  Distributed  Random  Variables:  A  wide 
variety  of  interesting  approaches  and  techniques  has  been  applied  to 
the  problem  of  finding  the  p.d.f.  of  a  product  of  n  normally  distrib¬ 
uted  random  variables.  In  fact,  as  a  collective  group,  these  investi¬ 
gations  represent  a  rather  extensive  study.  An  interest  in  this  problem, 
as  indicated  by  published  articles,  first  became  apparent  in  the  early 
1930's,  with  the  bulk  of  important  derivations  occurring  only  a  few 
years  later. 

In  certain  cases,  the  quadrature  method  used  in  treating  the  ex¬ 
ample  under  Theorem  3,  may  be  used  to  derive  an  expression  for 
Yf  =  N^,  the  product  of  two  normally  distributed  random  variables. 

These  special  cases  involve  the  dependent  and  independent  cases  of 

and  N2  described  by  the  normal  p.d.f. 'a  with  parameters  N(0,  1)  or 
2  2 

N(0,  )  and  N(0,  ).  These  four  results  are  enumerated  in  Table  I 

Appendix  A. 

The  dependent  and  independent  cases  of  Y^  =  in  which  N^, 

are  normally  distributed  with  arbitrary  means  ^  4  0)  and  variances 
do  not  readily  lend  themselves  to  the  ordinary  integral  methods, and  as 
a  result  have  undergone  extensive  investigation, 

A  chronological  history  of  the  important  results  in  the  study  of 
this  problem  is  outlined  in  order  to  suggest  certain  "difficulties"  in 
the  applications  of  these  results. 

In  1932,  Wishart  and  Bartlett  [l33]  considered  the  problem  of 


determining  the  p.d.f.  of  the  product  (N^  The  problem  was  posed  in 


the  following  framework. 


Let  X.  and  X-  be  two  normally  distributed  random  variables, 
1  d 

2  ? 

N^(0,  0^),  N^CO,  o^’)i  which  satisfy  the  joint  p.d.f.. 


^x1,x2)  = 


-  - j  exp - 

2  ^  °i  °2  V  2( 


2  2 

a _  /  V  2 !  Vg  ,  *;■  V 

i-p2)  Vo}2  °i  °2  ^2/i 


The  characteristic  function  of  the  product  X^  X 2  is  by  definition 

_  f  it  X-,  XJ 
E  [e  1  2J  ,  or 


f  r  -it  XX 
lX2(t)  85  /  /  e  ^(*1,  x2)  dxj  dx2 


(4-1) 


1 1-2  i  ? 


cr  t  +  (l-£  )  o: 


,2  rr-2  f2  ^ 
1  °2  ‘  | 


The  generating  function  of  the  semi-invariants  of  y^sX^^  :*-s 

K  =  ->i  log  |  1-it  0^  <5^  (  f  1)|  -Ka  log  1 1-it  cr2(  -  1)| 


(4-2) 


The  semi-invariants  were  so  named  by  Thiele  and  were  later  called 
cumulants  by  Fisher,  if  the  logarithm  of  the  characteristic  lunation 
(an  in  (4-2))  ,  a  function  of  t,  is  expanded  in  a  power  series  of  (it) 
which  converges  for  some  range  of  t  containing  the  origin  as  an  interior 
point,  the  somi-mvariants  k  are  the  coefficients  of 

v  /  v 

(it)  'vl  in  the  expansion.  A  simple  relation  exists  between  the  kv(y^) 

and  the  moments  of  Y^;  the  first  two  being  k^(y^)  =  E  *  k2(y^)  ”  ^ar  t^ll' 


e 


The  semi-invariunts  of  y^,  kv(y^) ,  are  the  coefficients  of  (it)v/Vl 


in  the  power  series  expansion  of  K.  Thus 

kv(yi)  =  Mv-l)i  a?  c^v  |  (p+Dv  t  (^-l)' 


(4-3) 


The  moments  of  the  product  may  be  investigated  through  the  appli¬ 
cation  of  (4-3). 

In  the  case  of  a  sum  of  independent  products  taken  from  the  joint 
p.d.f,  ^(x^,  x2^>  denoted  by  2=  +  •••  h  this 

relation, 


m 


K2  =  nK, 

holds  iso  that  from  (4-3),  the  semi-invariants  of  may  be  obtained 

from 

kv  =  ^  cr 'i  °2  nl_v  |(f  +  DV  +  (^  -  i)v|  • 


The  relation  of  (Jj.^  (t)  to  f(  2  )  satisfies 

.  (t)  =  ()*y  (t)j  =  |  1-2  i  £  o'£  o“  t  +  (l-(>2)  or2  a^J 


-yi  n 


= f  eltZ  f(  2  )  dS  . 
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Ey  inverting  this  characteristic  function,  it  follows  that 


f(S)  - 


2  TT 
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,  ,,  2.  2  2  . 2)720 
|l-2i  f  a~2  t  +  (1-p  )  <r2  t  ( 


dt. 
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(4-4) 
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becomes 


p  e  y±  iv  w 

(1-p  )  e  C  i  +  oo  iylw 

^<y  >  -  -  f*  - - ^-77-  dw 


2fi* 


/ 

p  i  -  oo 


c'nvty*” 


after  these  substitutions  are  made, 


(1“Q2)  <5^  <n;  t  =  £  i  -  w,  (1  -  <*;  <n,  =  X 


(4-5) 


From  the  theory  of  the  complex  variable,  it  is  clear  that  the 
contour  may  be  deformed  into  a  real  axis  as  j  £>  |  ■*  1J  thus  (4-5)  may 
be  written  as 


^(ll'  “  T 


(1-  )  e 


oo 


/ 


cos  y.^  w 
(l+w2/zn 


dw 


n  a,*4”  Pyi 

(1  -  f  )  e 


K"(n-1)  „  (v  ) 

yi  S(n-l)  tyl 


,  -oo  «=.  yx  oo 


=  **  2Un-ly j  yl 

where  ^(n  -jj(yi)  is  a  m°dified  Bessel  function  of  the  second  kind 

x  y  '  -  yl 


of  >Kn-l)J  order*  Fox1  n  =  or,  =  o"‘,  =  1  and  r  =  0  the  p.d.f.  of  yn  reduces 


to 

^(yx)  =  ft  Kq  ( ly-jl  )  ,  ya  /  0  .  (4-7) 

Wishart  and  Bartlett's  method  may  be  used  to  show  that  Z,  the  product 

X1  X2 

of  two  normal  variables  where  Z  =  —  «  r=r 

*1  ^2 


is  described  by 


which  possesses  a  singularity  at  z  =  0. 


Shortly  after  the  publication  of  this  result,  P.  T.  Yuan  [^3^] 
demonstrated  that  if  X^  and  are  independently  and  lognormally  distrib¬ 
uted,  the  product  Z  =  (X^  -  a)(X^  -  b) ,  where  a  and  b  are  the  upper 
(lower)  limits  of  the  range  of  X^  and  X^,  is  distributed  as  (^-7). 

The  analysis  of  G.  C.  Craig  [20],  [2l]  ,  following  in  1936 »  is  per¬ 
haps  the  most  notable  concerning  the  product  of  two  normally  distributed 
random  variables.  Craig  considered  the  bivariate  normal  p.d.f.  with 
parameters  (  )x  ,  m  ,  <T  ,  <T  and  coefficient  of  correlation,  p  ). 

X  —  Xl-l  X  n  \ 

X  d  1  d 

3y  Theorem  2,  this  joint  normal  p.d.f.,  under  the  transformation 
w  -  x1  x^s  y.^  =  y2t  may  be  written  as 


^<w,  y2)  = 


21T  er 
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In  an  effort  to  generalize  this  result,  Craig  introduced 


V1  "  ~ 

x. 


*i 


and 


y  „  — £ 

2  cr 
x. 


coefficients  of  variation,  and  considered  the  transformation  y..  = 


>°2 

,  the  reciprocals  of  the 

X1  X2 
1_tri  °2 


1 

x  =  ~  .  The  marginal  p.d.f.  of  y  ,  under  this  transformation,  may  be 
i  1 

expressed  us  the  difference  of  two  integrals,  Thus, 


^yl^  ~  Il^yl^  “  IP^yl^’ 


in  which 


Ijtyp  ■= 


2  A  2 
V1  “  2  P  V1V2  +  V2 

2(1^) 


2T 


('i-6) 
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4- 


(V  -  £  v2)  X  +  (v2  -  Q  V  ) 


dx 

x 


1 


and  I2(y^)  is  the  integral  of  the  same  function  defined  on  the  interval 
( -oo  ,  0) , 

In  an  effort  to  simplify  any  numerical  calculation,  Craig  reformulated 
(4-9)  as  an  infinite  series.  The  infinite  series  expansion  is  derived 
from  the  joint  p.d.f.  of  variables  U  and  Z  ,  where 


=  z. 


Under  this  transformation,  Equation  (4-8)  becomes 
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1-p 

?u)  =  Vrr  exp 
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The  term  exp 
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. 2  U 

‘ 

V 1-0 

\i-e. 

u 

in 

powers  of  u 

for  all  u,  u  / 

^  ’  ■ 

may  be  expanded  in  a 


(4-9) 


simplified  to  some  extent  by  substituting 


vev2  _  .  v2  -  e  v2  . 

a  R.,  and  —  , — =  K2  • 


In  the  expansion,  the  coefficient  of  ur  \  r  >  1,  is  (R-]R?  z)  ■ 


in  whioh^]r(>  •  ),  the  confluent  hypergeometric  function  of  order  r,  is 

^2  T-,  _^3 

(r+2)vc'2!  <r+3)w/3» 


*1  olO 


R_  R_  z  (R1  R~  z)2  (R-,  R?  z)^ 

X  (R.  R_  z)  =  1  +  ■  ’  /by-  +  - m  +  •  •  *  , 

r  12  r+1  oi  f  ^  ^  X  ! 


with  (r+k)^  s  (r+k)(r+k-l)  .  .  .  (r+1)  . 

2 

By  this  expansion,  the  p.d.f.  of  Z  -  /  <3^  <31,(l-p  )  may  be 

expanded  in  an  Infinite  series  involving  confluent  hypergeometric  functions 


^  These  functions  are  discussed  in  detail  in  Whitaker,  E.  T.,  and 
G.  N.  Watson,  A  Course  in  Modern  Analysis,  Cambridge  University  Press, 
Cambridge,  195^ 
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powers  of  z,  v^and  This  series  is 


*,  ,  i/I  -  P2  i  ''I2'2  ?  Vs  **2  e z 

^(z)  exp  < - p - +  - 5~ 

^  (  2(1-5  )  U-5  ) 


1 2; 


2T2 


If  (RjR,  z)  K1(z)+(H1\r24)  ^  (R^  z)  K2(z)  + 


+  (R16+R26)  'flrZ^6  (H1R2  z)  K3(z)  +  *  *  *  ♦ 


where:  K^(z)  =  the  Bessel  function  of  the  second  kind  of  the  ilb  order 
and  with  argument  z,  and 


Ep1*'’-#  (i>  ^VvT-’*- 


which 


I^(  •  )  is  the  Bessel  function  of  the  first  kind  of  the 
jib  order. 

When  v1  =  v0  =  o,the  p.d.f.  of  Z  =  X2  /  is  the  simple  Bessel 

function  result  expressed  by  (4-7). 

Craig's  results  have  unfortunately  proved  to  be  of  little  use  com¬ 
putationally  for  it  may  be  shown  that*  for  large  v^  and  v,,  the  series 
expansion  (4-9)  converges  very  slowly;  in  fact  for  v^  and  v2  as  small 
as  2,  the  expansion  is  unwieldy.  Yet  after  this  publication,  many  in¬ 
vestigators  consider  the  problem  to  be  "solved"  although  the  convergence 
problem  exists.  Others  have  been  critical  of  this  "solution",  even 


Craig  himself  who  in  19^2  ,  [2lJ  ,  stated  that  even  thougn  his  result 
is  a  mathematical  solution  to  the  problem,  it  falls  far  short  of  what 
is  required  for  numerical  computation. 

Using  Craig's  formula  for  the  cumulant  generating  function, 

J.B.S.  Haldane  in  19^2  presented  the  moments  (about  the  mean  and  the 
origin)  and  the  cumulants  of  various  products  and  powers  of  normally 
distributed  random  variables.  He  considered  the  correlated  and  un¬ 
correlated  cases  of  these  statistics:  a)  the  cube  of  a  normal  variable, 
b)  an  arbitrary  power  of  a  normal  variate,  and  c)  the  product  of  n  normal 
variates.  These  expressions  are  rather  cumbersome  however  and  their 
application  in  an  applied  problem  requiring  an  extensive  amount  o.f  numerical 
computation  would  be  severely  limited.  This  fact  combined  with  the 
realization  that  Craig's  series  is  subject  to  restrictive  conditions 
prompted  others  to  examine  this  product  more  from  the  viewpoint  of  establishing 
its  analytical  properties.  For  instance,  B.  Epstein  [31 ]  and  I.  Kotlarski  [&8j 
independently  showed  that  theMellin  transform  is  a  very  useful  analytical 
tool  in  examining  the  integral  equation 


'l>  ■  j  \  \  ^2  >  d?2  • 


Their  investigation  are  closely  related  to  problems  of  the  nature  of 


those  discussed  in  Chapter  VI. 


In  19^7 ,  L.  Aroian  [2]  took  up  the  problem  of  convergence  in 
Craig’s  series  expansion.  Using  Craig’s  notation,  he  showed  that  as 


and  V£ - *■  00,  the  p.d.f.  of  Y^  approaches  the  normal  p.d.f.  In 

addition,  he  demonstrated  that  the  Type  III  function  and  the 
Gram-Charlier  type  A  series  afford  excellent  approximations  to 
the  distribution  of  Y^  when  ^  =  0. 

Using  the  properties  of  the  moment  generating  function  of  the 
^(y^) ,  it  is  possible  to  show  that  £i  [y^]  =  Yp  =  Vg  +  6  and 


the  standard  deviation  is  <T 

y 
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Aroian  proved  the  following  statements  in  the  form  of  theorems: 


1)  The  p.d.f.  of  Y^  approaches  the  normal  p.d.f..  with  mean  y^ 


and 


and  variance  <r  as  V,  and  V,- 
y,  1  2 


00 (  or  -00 )  in  any  manner 


whatsoever,  provided  -1  +  €  *  1,  £  >  0. 

2)  The  p.d.f.  of  y^  approaches  the  normal  p.d.f.  with  mean  y.^ 

2 

and  variance  if  - *•  00,  V - *»  -00,  provided 

-1  *  ?  *  1  -  €  ,  €  0 


3)  The  p.d.f.  of  y^  approaches  N(y  ,<r^  )  if  remains  constant  and 

1  yl 

- *r  oo,  -1  +  e  ^  R  —  1,  €.  >  0;  or  if  remains  constant 


and  V. 


-00  for  -1  ^  R  <  1  -  e  , €  >  0. 


Aroian  demonstrated  the  close  approximations  to  P^(v1)  by 


5.5 


W>j 


PW> 


the  Type  III  function  and  the  Gram-Charlier  type  A  series  by  numerically 
integrating  Craig's  expression  (4-8)  for  the  special  easel =  0, 

=  10,  and  p  =  0.  A  brief  tabular  comparison  is  presented  in  hia 
article. 

In  Table  X.  ,  Appendix  A,  the  numerical  integration  of  Craig's 
formula  (4-8)  for  a  few  special  cases  is  presented These  are; 

(1)  Vx  -  V2  »  f  -  0V  (2)  Vx  =>  V2  =  %  ?  =  0  and  (3)  V;L  «  1(0), 

v2  =  0(1),  p  =  0. 

4.6.2  The  Quotient  of  Normally  Distributed  Random  Variables:  The 
first  investigations  of  the  properties  of  the  p.d.f.  of  a  quotient  of 
two  normally  distributed  random  variables  were  directed  toward 
characterizing  the  quotient's  properties  in  terms  of  the  properties  of 
the  component  variables. 

K.  Pearson's  study  [97]  in  1910  of  an  opaonlo  index  formed  by 
the  quotient  of  two  normally  distributed  random  variables  represents 
the  first  published  investigation  of  this  problem.  He  succeeded  in 
obtaining  the  first  four  moments  of  =  X^/x^  in  terms  of  the  moments 
of  X^  and  X Unfortunately,  he  found  that  they  were  "practically 
unworkable  if  X^  and  X2  are  correlated  as  we  should  have  to  find  the 
third  and  fourth  order  product  momenta". 

Later,  C.  C.  Craig  [20]  ,  (1929) *  developed  this  approach  by 

IT) - 

A  number  of  other  cases  are  to  be  included  in  a  forthcoming  research 

report  to  be  released  ly  the  Applied  Mathematics  Research  Laboratory, 

Aerospace  Research  Laboratories. 
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deriving  the  moments  of  Y^  in  terms  of  the  semi-invariants  of  the  eom- 
ponents.  Craig  demonstrated  the  advantages  of  using  the  semi -invariants 
(as  opposed  to  the  moments)  by  constructing  the  momenta  of  Y^  in  the 
oase  of  correlated  ^  and  with  relatively  little  difficulty.  The 
resulting  expressions  are  cumbersome  by  present  day  standards  and  are 
difficult  to  apply. 

In  this  time  period,  the  hypothesis  that  Y.^  must  be  near-normally 

distributed  given  that  the  components  X.^  and  X^  are  distributed 

normally  was  strongly  supported.  Investigations  of  Merrill  [86]  in 

1928  aud  by  Geary  [38]  ,  1930,  were  instrumental  in  disproving  this 

conjecture.  Merrill's  investigation,  by  graphical  approximation, 

showed  that  when  the  correlation  between  X,  and  X„  is  high  and  the 

or  j.  2 


coefficients  of  variation  ^x.  =  -  * 

i  ji 


,  i  =  1*  2,  are  large,  there  1b 


a  considerable  deviation  from  the  normal  p.d.f.  Geary  established 
this  result  on  a  more  rigorous  foundation  by  formulating  what  is  now 
a  widely  known  approximation  He  considered  the  problem  as  formulated 
in  this  manner: 

Let  X^  and  be  two  jointly  distributed  normal  random  variables 
with  the  p.d.f., 


g(x1,x2)  = 


jt  txp  [  (c?)  *&L^n**i 


(4-11) 


where : 


~  Pi'  ~  P2'  =  =  G~2.  aQd  ?  is  the  of 
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correlation 


When  the  expected  values  =  0,  the  gC*^*  reduces  to 


g(xltx2) 


2^01 


re  H 


exp 


2(1-^) 


2  * 


°i 


X1X2  +  *2  \ 

°2  °i  / 


.  (4-12) 


Consider  the  random  variable  Y^,  where 


*1  = 


b  +  x. 


a  +  x. 


(4-13) 


and  a  and  b  are  constants;  then  the  function 


a  y  -  b 

t  a  1^===-::^;=  (4-14) 

V  yi2 ' 2  «  °i  °2  h  * 


is  approximately  normally  distributed  with  expected  value  t  =  0  and  standard 

deviation  <f£  =  1,  provided  that  (a  +  x^)  is  unlikely  to  assume  negative 

values.  Geary  shows  that  this  latter  condition  is  satisfied  if 

X 

(a  »  3  <J^),  i.e.,  the  coefficient  of  variation  of  (a  +  x^)^j  . 

Geary's  and  herrlll's  results  gave  definite  proof  that  tho 
hypothesis  of  the  "normally  distributed  index"  was  incorrect  in  many 
oases  and  have  lead  to  the  problem  termed  ''characterizing  properties  of 
a  p.d.f."  which  is  discussed  in  Chapter  VI. 

It  is  advantageous  at  this  poiht  to  demonstrate  that  the  integral 
methods  dlsoussed  in  Chapter  III  may  be  used  to  obtain  an  expression  for 
the  p.d.f.  of  in  order  to  lay  further  groundwork  for  discussing 
other  known  results. 


In  the  independent  case*  the  joint  p.d.f.  of  and  is  of  the  form 

g(xlf  x2)  *  f^x^)  f2^x2^*  C 4-15) 

When  considering  the  random  variable  Y^  =  X^/X^,  it  is  necessary  to  use 

the  methods  described  in  Section  3 In  order  to  define  a  one-to-one 
transformation.  By  specially  defining  Pr(X2  =  0)  =  0,  the  inverse 

functions  x^  =  y^/y2  and  x2  =  y2  map  the  set 

V  e  j  -oo  <  x^  *  oo;  -co  c  x2  ^  0,  0  «£  x2  <  oo  \  into  the  uet 

A  -  {  -oo  ^  y1  *  oo  ;  -oo  -*•  y2  0,  0  <  y2  ^  oo j  where  the  joint  p.d.f. 

*<yr  y2)>  0.  The  I  J  I  =1  y2 1  and  the  marginal  p.d.f.  of  Y^  is,  by 

equation  (2-27), 


Equation  (4-l6)  is  directly  integrable,  provided  the  limits  exist,  so 
that 


g(y  )  =  — i  - - ,  y  ^  0.  (4-17) 

ir(yi2  <r22  t  of) 
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The  more  complicated  result,  /  u2  /  0  of  g(y^)  may  also  be  ob¬ 
tained  by  applying  the  methode  of  Seotion  3 . 4.2.  This  result  was  first 
shown  by  Heller  [ 33 J  and  Baker  [4  3  in  1932. 12  It  is 


gCy,)  “ 


ri  °2  f' 


1  ’  W  -i 


P-l  -2  Wz  »2 

+  — 

°i  °2 


(^1  -  *2y 


ai(*J*i  °i>+yi  <r2) 


°22  *i2  - 2 ^ +  °i2/f '  ^"a2  yi2-2  **i  °i  °i +  3/2 


/ 


°g  ~  +  y-^  )ip  <3^  -  jUi  <r2) 

~1  °2  <^><42  *7  -  2?*X  «i  “5  +  'ri2)  * 


iA  2 

e"A  u  du  (4-18) 


We  see  that  if  ^  s  ^  s  ^  s  0,  g(y.^)  reduces  to  Eqxiation  (4-17);  and 
under  the  conditions  =  0,  p  ^  0,  g(y.j)  reduces  to 


«(y1) 


yi-?  "l  ■’I _ 

^Ce?2  y2  -  2  ^  crx  <r2  yx  +  <^2)  ** 


(4-19) 


The  calculation  of  values  of  the  c.d.f,  of  y1  is  accomplished  merely  by 


ip) 

This  result  is  generally  attributed  to  E.  C.  Fieller. 
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the  use  of  a  set  of  tables  of  arc-tungent  in  these  two  oases. 


Fieller's  paper  is  mainly  devoted  to  obtaining,  by  quadrature,  an 
expression  which  will  simplify  the  numerical  calculations  associated 
with  equation  C 4— 18) .  His  rather  complicated  quadraturo  method, 
formulated  in  terms  of  existing  tabulated  functions  of  the  bivariate 
normal  density  function,  was  derived  under  the  following  hypothesis. 
Consider  the  ratio  y^  =  x^/x2  where  random  variables  X1  and  X2 

are  described  by  the  joint  p.d.f.  of  equation  (4-11).  The  points 
(*1’  x2^  ^responding  to  a  given  value  of  lie  on  the  line 


*1  =  y-^  x-  .  (4-20) 

Thus,  the  probability  that  an  arbitrary  element  of  g(x^,x^)  will 
have  an  index  {  vj  *  <  V2  [  equal  to  the  volume  of  the  portion  of 

g(xL,  x2)  which  lies  above  the  area  swept  out  in  the  (x^,  x^)  plane 

by  the  line  x1  =  y^  x^  as  it  revolves  from 


When  g(x^  ,  x^)  is  of  the  fo"m  (4-ll)  and  each  variable  is  expressed 
in  terms  of  deviations  from  its  respective  mean,  V  »ay  be  written  as' 


v  =  f  f  -  1  ! 

(  1 

r  2  ,  2  1 

x-  2rx,x_  x, 

2  12  t  1  f 

J  j  2rr<r,<r2V^2  j 

a+b  1  * 

1  2(l-r2) 

^  2  2  ( 
^2  °1  °2  °1  ' 

dx1dx2 


(4-22) 


where  a  and  b  are  two  portions  of  the  (x^  ,  x^)  plane  indicated  in 
Figure  3«  The  boundaries  of  the  a  and  b  portions  are  the  lines 

x2  +  x2  " 

*1  +  *1  »  yl^x2  +  X2^* 

If  the  transformations  =  o(  and  x^  -  y^  x^  *  71  are  made,  the 

portions  X  and^Q  of  the  (  o(  ,T]  )  plane  which  correspond  to  a  and  b 
in  Figure  3  are  bounded  by  the  lines 

o<+  x  =  0, 

7|+  ^  -  yx  x2  =  o. 

Th±6  area  is  shown  in  Figure  4, 

Therefore  if  this  change  of  variable  is  performed,  V  becomes 


-f  f _ i— 

j  j  ?TT  VII 


-y?x  tr- 

«  "  do<d  7| 


where  x 


2r  (o0(7]+  y^  o<  ) 


_1 _ (  _j£_  2r  (c0(T 


(  ^  +  y1®( ) 


(4-23) 

,.2  \ 


n- 


In  this  derivation  r  denotes  the  coefficient  of  correlation. 


Figure  4.  Geometrical  Presentation  of  The  Bounded 
Area  in  the  Transformed  Case. 
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which  is  identically  - =- 

X-P 


-°L  XL 

<r  <r 
oL  \ 


From  this  identity,  these  results  are  apparent s 


*£  (l-e2)  =  4J2  er/  d-r2)  -  2  py1  <JJ  +  y/  dg2) 


(l-^2)  =  ffj2  (1-r2) 


(4-24) 


*2  =  (r  dj  -  y;L  dg)  /  og  <r^~  (1-r2)  . 


By  squaring  the  last  equation  of (4-24)  and  multiplying  the  result  by 
the  first  two  equations,  we  obtain  these  results: 

£2  ~  ”  yl  °z ^  j  ^dl2  ~  2  r  yl  °i  ^2  +  yl2  tf22'>’  (^-25) 


so  that 


(l-(>2)  =  (1-r2)  of  /  (of  -  2  r  61  (Sg  +  y^  of),  (4-26) 


\’ei> 

^  -  (d^2  ~  2  r  yl  ^1  °2  +  °22^  ’ 

<T ^  <5^  ^fl^2  =  <5^  <5^  y  1-r2  . 


If  we  write  X  = 


o4  =  ^2 

V  *5 


(4-27) 


(4-28) 


(4-29) 
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and  1  = 


m 


$ 

IS 


$1 


$ 


:s 


I 


*  *1  "  yi  x2 


2  _2* 


>(  <«i  -  2  r  yi  °i  °i  +  yi  ^  > 

the  quadrants  A  and  B  of  the  (X,  Y)  plane  that  correspond  to  the 

portions  a  and  b  of  the  (x^,  x^)  plane  have  as  a  common  corner  the 

point  (-h,  -k)  where 


*2 

h  =  — , 


.  X1  ~  yl  *2 _ . 

(<Ji2  - 2  r  yi  °i  ^  +  yi2  *22) 


(4-30) 
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From  equation  (4-23)»  V  becomes 


=  f  f  ■  - - -  exp  -I-  r  (X2  -  2  o  XT  +  f2)  dX  dY  L  (4-31) 

J  J  ZlTcrcr^J  1-r  5  2(l-p?) 


the  probability  of  obtaining  an  index  for  which  y^  3s.  V  is  given 


oo  oo  oo  oo 


1  -  V  =  /*  /*  +  f  f  - jL-k  exp  “  2  fx2  -  2  e  XY  +  Y^dX  dY . 

/  *  &  J*  2ti-p  > 


(4-32) 


Equation  (4-25)  provides  two  values  of  The  appropriate  value  is 
chosen  by  noting  that  as  y^  - 1>  oo  ,  the  point 


(h  ,k) 


* 


* 

i 

i 


-  x. 


so  that,  the  positive  square  root  makes  (1-V)  - — >  0  as  is  desired. 
Extensive  tables  are  found  in  [l4o]  for  whioh  the  value  of  the 
integral, 


h  k  h  k 


is  calculated  for  small  increments  of  ^  ,  -1*^>  41.  For  non- 
positive  h  and  k,  the  following  relations  hold: 


I  (2) 

I 

I 

and 


ffzW 

h  -k 


dXdY 


//  Z(-^)  dXdY 
h  k 


Probabilities  of  interest  may  be  calculated  with  the  aid  of  these 
relations  and  the  appropriate  values  of  p,  h,  k,  X,and  Y. 

■ith  the  publication  of  Heller' b  quadrature  method,  many  agreed  that  the 
problem  of  the  quotient  of  two  normally  distributed  random  variables 
was  essentially  solved.  In  later  papers,  emphasis  was  turned  to  the 
study  of  the  general,  mathematical  properties  of  the  random  variable 
quotient.  Two  such  notable  papers  appearing  in  the  late  1940 'a  were 
published  by  Curtiss  [24]  and  Gurland  [48]  .  Curtiss  approached  the 
study  of  the  properties  of  the  quotient  from  an  analysis  of  the 
application  of  the  Lebesque-Stielt jes  integral. 

Gurland  formulated  the  c.d.f.  of  the  ratio  of  two  linear  combinations 
of  correlated  normal  random  variables.  He  presents  two  theorems  for 
this  distribution: 

Theorem  1:  (for  positive  denominator).  Let  X2,  .  .  .  ,  Xq  have  a 
joint  c.d.f.,  F(xx ,  r  ^  •  •  *  ,  xfl)  the  corresponding  characteristic 

function  t 2,  .  .  .  ,  t^).  Let  G(x)  be  the  c.d.f.  of  the  linear 

combination,  (a.,  X,  +  a„  X,  +  .  .  .  a  X  )  /  (b.  X.  +  ,  .  .  b  X  ), 

11  d  z  nn  ii  nn 

where  a^  a2>  .  .  .  ,  aR,  b^  ...  ,  bQ  are  real  numbers. 


Theorem  2:  (for  positive  or  negative  demoninator) .  Let  G(x)  be  the 
c.d.f.  of  the  ratio  of  linear  combinations  (of  Theorem  1).  If 


Ebi*-°|-o 

i=l  ' 


a(x)  +  g(x-o) 


/j*  it(vbi  ^ . ‘’vWl  ♦ 

t 


j-bl  X11 . t(Vbn  Vi 


where:  /+(t1,t2 . tn)  =  ///•‘ 


(4-35) 


( t-  x. f < • » |t  x  )  ,  v 

11  n  n  d  F(x, ,x_, . . . ,x  ) 

12’  n 


E",.  *k>0 


0  <tlft2,...,tn)  =  ff'j 


i(t  x  ...  t  x  ) 

0  d  F(x^,x2 * . . . ,xq) . 


y>  X.  ^0 

'  n  k- 


4.3  The  Rectangular  Density  Function 

A  random  variable,  X,  defined  on  a  finite  interval  =  {  a  <  x  6b  |  , 
(a,  b^ finite  real  numbers) f is  said  to  obey  the  uniform  or  rectangular 
p.d.f.  if  over  the  finite  interval  "if  ,  the  probability  of  A,  a  sub¬ 
interval  of  if  ,  is  given  by 


ifSHH  * k  a  subset  °f  r 


O  ,  otherwise. 


(4-36) 


i'lTcl'ii 


mm 


Thus  the  c.d.f.  of  the  rectangular  probability  law,  by  this  definition,  is 

10,  x  -«  a, 

~  ,  a  =  x  2  b,  (4-37) 

1,  b  •*  x. 

If  we  differentiate  (4-37),  we  find  the  p.d.f.  of  the  rectangular  error 
function  to  be 


f(x)  = 


b-a 


=  0, 


x  2s  b, 


elsewhere. 


(4-38) 


4,3.1.  The  Product  of  Two  Identically  Distributed  Rectangular  Variables; 
The  product  of  two  identically  distributed  rectangular  variables,  Bay 

y^  -  R.^,  easi1Jr  found  by  quadrature  utilizing  the  methods  of 
sections  3.4.1.  and  3-4.2.  Four  distinct  oases  arise  in  solving  for 
the  p.d.f.  of  =■  R.^,  wliere  ^  a  rectangular  vuriable  defined  on 
US  [a  ^  x  ^  b  |  .  The  distinct  cases  may  be  enumerated  as:  Case  (1) 
y^  =  R^R^,  ^  defined  on  interval  ^  so  that  0  ««  a  ■**  b,  Case  (2)  R  de¬ 
fined  on  interval  &  so  that  a  0  b,  |bl  >•  |  a  |,  Case  (3)  R  de¬ 
fined  on  interval  so  that  a  -*■  0  b,  |aj  s-]b|i  and  Case  (4)  R 

defined  on  interval  V  so  that  a  -*•  b  -*  0.  The  p.d.f. 's  of  y^  for 
each  of  these  four  caseB  are  presented  in  Table  II,  Appendix  A.  The 
distribution  of  the  product  of  any  number  of  identically  distributed 
rectangular  variables  can  easily  be  obtained  from  Kendall's  jl46 
derivation  of  the  p.d.f.  of  the  geometric  mean  of  samples  from  a 
rectangular  distribution. 

4.3.2,,  The  Quotient  of  Two  Identically  Distributed  Rectangular  Variables: 


The  quotient  of  two  identically  distributed  rectangular  variables  is 


also  easily  found  using  the  methods  of  Sections  3.4.1  and  3.4.2.  The 
same  distinct  cases  depending  upon  the  Interval  of  definition  as  in 
4.3.1  arise  in  the  problem.  The  p.d.f. 's  of  theso  four  oases  are  pre¬ 
sented  in  Table  III,  Appendix  A. 

4.4  The  Triangular  Density  Function. 

The  sum  or  the  difference  of  two  identically  distributed  reotangular 
defined  on  |  a  *  x  £  b  j  is  triangularly  distributed.  The  p.d.f. 
of  T  is 


f(x) 


Ix-2  a  _  , 

- 2  a  <  x  <■  b+a  , 

(b-a)2 

— — b+a  <  x  <  2  b  , 
(b-a)2 

0  ,  elsewhere. 


(4-39) 


4.4.1  The  Product  of  Two  Identically  Distributed  Triangular  Variables: 
Slightly  more  cumbersome  expressions  result  from  applying  the  methods 
of  Chapter  HI  to  finding  the  p.d.f.  of  Y^T^  and  .  This 

difficulty  is  due  primarily  to  the  fact  that  the  p.d.f.  of  the  tri¬ 
angular  variable  is  defined  by  two  branches  as  in  (4-39).  In  deriving 
the  p.d.f.  of  T^,  one  must  determine  which  combination  of  branches 
of  Tj,  and  must  be  considered  for  each  of  the  various  mutually  ex¬ 
clusive  partitions  of  the  set  defined  under  the  transformation.  The 
condition  is  best  illustrated  by  considering  the  problem  of  deriving 
this  p.d.f.  when  T  is  defined  on  an  arbitrary  positive  interval.  Under 
the  transformation  y1  =  ^  =  x1  x^  ,  y0  =  x,  ,  the  set 


V  =  |  2a  -*  -«  2b,  2a  x^  -«  2b  j  is  mapped  into  the  set  graphically 

illustrated  in  Figure  5.  Here,  the  notation  g,  g*  and  h  ,  h'  denote  the 
respective  branches  of  f(x^)  and  f(xg)  which  apply  to  the  partitions  of  (i  . 
The  p.d.f.  of  is  derived  by  integrating  over  the  product  combinations 
as  indicated  in  Figure  5  over  the  partitions  A,  B,  C,  D,  and  E  of  /3  . 

Two  important  special  cases  have  been  derived-  (Table  IV,  Appendix  A). 

These  are  the  cases  in  which  T  is  defined  over  a  positive  interval  ^0,  2  j 
and  t'.e  interval  J^-2,  2~j  . 

4.4.2  The  Quotient  of  Two  Identically  Distributed  Triangular  Variables : 

The  p.d.f. ' s  of  two  opecial  cases  for  T  defined  on  ^0,  2j  and  ^-2,  2J 
have  been  derived  and  are  presented  in  Table  V,  Appendix  A. 

4.5  The  Product  and  Quotient  of  "Mixed  Components11. 

The  product  and  quotient  forms  of  the  three  measurement  error 
variables  under  consideration  which  have  not  yet  been  treated  are 
RN  ,  TN  ,  R/N  and  ‘f/N.  There  is  no  requirement  to  investigate  the 
reciprocal  forms  of  these  quotients  since  Cramer  (^22  j  has  insured 
that  the  p.d.f.  of  the  reciprocal  of  a  quotient  is  immediately  evident 
upon  derivation  of  the  p.d.f.  of  the  quotient. 

The  four  forma  above  lead  to  rather  cumbersome  expressions  in  re¬ 
spect  to  numerical  calculations.  A  more  convenient  method  of  analysis 
in  the  cases  of  the  quotient  is  demonstrated. 
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4 *2  2a(b+a)  4ab  (b+a) 


2b(b+a) 4b2 yl 


Figure  5*  Geometric  Representation  of  the  Region , 


4.5.1. 1  The  Product  of  a  Reotangular  and  a  Normally  Distributed  Variable: 

Certain  cases  of  this  random  variable  product  may  be  formulated  in  terms 

of  the  tabulated  exponential  integral.  As  an  example,  consider 

=<  RN  =  Xg  where  X^  is  rectangularly  distributed  on  (0<x^  SL  1  j 

and  X2  is  normally  distributed  with  parameters  N(0,  1)  and  defined 

such  that  PrCX^  =  0)  =  0.  Defining  the  transformation,  y^  =  x^, 

y2  =  ,  the  p.d.f.  of  Y^  is  by  Section  3.4.2,  Chapter- III, 

2 


f  *x  ,  .  ■* 
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-00 


1  e_ 


dy 


2  ’ 


-00 1  yn  £  0, 


^(yx)  =< 


oo 

/ 


1  e 
yjTn 


-fe  y5 


dy2  » 


9  <  00 . 


V. 


Thus 


,  with  the  substitution.  <T  =  yg  ,  ^(y^)  satisfies 


^(yJ  =  - 7=  / 

2  V  2TT  ^5 


00  -jfcr- 

e 


dcr- 


yx  ^  0 ; 


which,  in  turn,  may  be  expressed  us  a  function  of  the  tabulated 
exponential  integral  after  the  substitution 


v  =  }tcr  . 


s 


All  product  forme  of  Y^  =  RN  ,  where  N  is  taken  as  N(0,  1),  may  be 
expressed  as  a  function  of  the  tabulated  exponential  Integral.  An  extensive 

table  of  this  integral  is  presented  in  reference  [l 43  ]  . 

2 

Other  product  forms  in  which  N  is  N(yU  /  0,  <T/0)  are  not  ex¬ 
pressible  as  simple  functions  of  this  integral, 

4.5. 1.2  The  Product  of  a  Triangular  and  a  Normally  Distributed  Variable: 


As  was  the  case  in  the  preceding  section,  only  the  product  cases  of 
Y^  =  TN  in  which  N  is  taken  N(0,  1)  are  easily  attainable  by  the 
quadrature  methods.  One  of  the  less  complicated  examples  is  the  case 
of  Yx  =  TN  a  X±  X2  where  T  is  defined  on  [0  <  xx  <  ?.  j  and  N  ia  N(0,  1) 
The  p.d.f.  of  is, 


f(xl}  = 


g^)  =  x1  , 


g'(xx)  =  2-xx, 


0 


1  <  Xx  i  2  . 


The  p.d.f.  of  X2  is  defined  such  that 


f(x2)  - 


1  X2 

hCxJ  =  — ==.  e 

^  nr 


-oo  *  x2  <  0,  0  <  x2  <  oo 


-  0  , 


x2  =  0. 


The  regions  of  definition  under  the  transformation  y1  =  x^  x2,  ~  x2  1 

shown  in  Figure  6. 


4.5.2. 1  The  Quotient  of  a  Rectangular  and  a  Normally  Distributed  Variable: 
A  rectangular  variable,  say  Xq,  defined  on  |  a  £  *Q  b  |  be  re¬ 
defined  on  the  interval  [0,  l)  through  the  use  of  the  transformation, 

X  =  X(b-a)  +  a 
o 

where  X  is  rectangularly  distributed  on  [0,  l)  .  The  original  rectangular 
variable  may  also  be  redefined  on  the  interval  (-1,  l]  by  the  trans¬ 
formation 

*„  - 

where  =  half-range  of  the  interval  [a,  b]  and  >£  =  mean  of  the 
interval  [a,  b]  . 

The  .distribution  of  Y  =  R/N  is  greatly  simplified  by  raing  one 
of  the  above  transformations.  In  the  special  case  of  =  R/tt  =  Xj/Xg 
where  R  is  defined  on  {  0  £  x^  £  1  j  or  [  -1  £.l}  and  N  is  taken 

to  be  N(0,  1)  where  Pr(X2  =  0)  =  0,  the  p.d.f.  of  y1  is  simply. 


In  the  more  general  cases,  R  may  be  distributed  as  any  of  the  four 

2 

cases  enumerated  in  Section  4.5,  and  N  taken  to  be  N(  tf~  ).  The 
respective  p.d.f.' a  of  y1  in  these  four  cases  are  shown  in  Table  VII, 
Appendix  A.  It  is  apparent  that  these  expressions  are  rather  awkward 
and  difficult  for  numerical  computations. 

Broadbent  [15]  has  contributed  a  general  method  with  which  to  in¬ 
vestigate  a  ratio  in  which  the  numerator  is  a  rectangularly  distributed 
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random  variable  and  the  denominator  is  unspecified  but  possesses  certain 
general  properties.  The  c.d.f.  of  =  R/lN  may  be  obtained  in  terms 
of  tabulated  integrals  by  treating  this  special  case  by  Broadbent’s 
method  of  analysis. 

Any  arbitrary  random  variable  quotient  of  independent  components 
whose  numerator  is  reotangularly  distributed  may  be  written  in  the 
following  "standardized"  form: 


y^  =  (  V  )/(  /i  +  d*  £  ) , 


(4-43) 


where : 


Y[  is  rectangularly  distributed  on  [-1,  l]  , 

V  is  the  expected  value  of  the  original  R  variable, 
is  the  half-range  of  the  original  R  variable, 

§  is  normally  distributed  with  parameters  (O,  1), 

and  jx  and  6~  are  the  parameters  of  the  original  N  variable. 

The  hypothesis  is  only  slightly  restricted  by  requiring  that 

^  0,  _/S/  '€  =  oC  where  0  <  oC  <  1,  *  0,  —  -  ^  where  O*  j£  *  oo  . 


Since  y'  =  (  ™—  )  (1  +  taC  \  )/(l  +  ^  £  ) , 

X 

the  standardized  quotient  y^  may  be  considered,  where 


(4-44) 


(4-45) 


In  (4-45),  is  rectangularly  distributed  on  (1  -  o£,  1  +  oO  and 
ia  normally  distributed  with  mean  1  and  variance  .  In  this 


special  case,  *2  is  described  by  a  probability  density  function 

j 

f(x2)  =  —  F(x2)  existing  almost  everywhere  and  a  c.d.f.  F(x2)  =  z 
which  has  au  inverse,  say  G(z),  which  is  defined,  non-zero, and  finite 
for  almost  all  a.  This  condition  allows  the  preclusion  of  probability 
measure  at  points  (x2  =  0)  and  (x2  =  oo  ). 

The  function  z  is  independent  of  X^  and  is  rectangularly  distrib¬ 
uted  on  (0,1),  a  property  of  all  c.d.f. 's.  The  joint  p,d.f.  of  z 
and  is  by  theorem  2,  Chapter  ill, 

*•]_)  =  dz  dxi* 

In  order  to  obtain  the  joint  p.d.f.  of  z  end  x^,  note  that  for  any  z, 
dx^  =  |G(z)|dy^  and  y1  is  bounded  by 

((1  -  <*0/G(z),  (1  +  o6)/G(z)),  if  z  >C» 

((1  +  od)/G(z),  (1  “  oO/GCa)),  ifz-b-C, 

where  F(0)  =  C  ,  £  not  necessarily  small. 

Therefore,  z  and  y^  have  the  joint  p.d.f., 

yx)  =  §-^|G(z)j  dz  dy» 

over  the  region  shown  in  Figure  7> 


The  Pr(0  *  y^  £.  Q^),  for  (Q^  >0),  Ls  given  by  the  integral  of  this 

function  over  the  region  A  bounded  by  z  ~  1  and  y^  =  (1  -  a0/G(z) 

minus  the  integral  over  the  region  S  by  y^  =  Q^,  z  =  1  and 
y^  =  (1  +  e(.)/G(z),  so  that 

Pr(0  <  y;L  *  :=  f  f  ~~r  | GC a) |  dz  dy1  -  f  J  — — ■  |g(z)|  dx  dy^  (4-46) 

A  S  ^ 

After  integration  of  this  repeated  integral  in  respect  to  y^,  the  first 
integral  nay  be  expressed  as 


A  /  r  ,  K  «*>  -  (1  -  *>  }> 

z  =  F  [(1-oO/qJ 

If  we  change  the  variable  to  £  end  write 

g(f )  =  j*fa  i ), 


(4-47) 


the  result  is  the  p.d.f.  of  the  standardized  form  of  X^.  Now  let 
X  =  (1  -  <  -  Q-, )/(  V'Q.J, 


The  integral  (4-4?)  becomes 


V^Q, 

— -  J,  g(X) 


r 

J  X)  s(f  )  df  . 

r= * 


In  this  special  case  where  random  variable  X,  is  normally  distributed. 


1  i^H  >|  Ij _ 


g(X)  *  In(x)  ifl  Hermitian  probability  integral  of  order  n 
tabulated  in  the  British  Association  tables,  2nd  ed.  (1946)P 

The  second  integral  of  (4-46)  is  obtained  in  the  same  manner  and  is 
denoted  by  g(Y).  Since  Pr(y1  -£  0)  =  Pr(x2  ^  °)  +  Pr(x2  =  oo)  =£ 

as  previously  noted,  the  c.d.f.  of  yx  is,  for  Q  >  0, 

=  Pr(yi<  Q1)  *C  *  z-g  {j1  ft(X)  -  ^  g(Y)|  *  (4-48) 

It  follows  that  from  (4-48),  the 


Pr(y1  >  Qx) 


g(Y)]  -  [: 


X  +  J, 


(4-49) 


The  forms  (4-48)  and  (4-49)  ere  useful  when  £  is  small  for  the 
calculation  of  percentage  points  at  the  tails  of  the  distribution  of 

yi since  as  Qx  — »  0  and  X,  Y  - >  00 ,  both  ^  g(X)  and  J  g(Y)— ►  0  - 

the  former  more  slowly  than  the  latter.  For  small  Q1,  J  g(Y)  is  small 
in  comparison  with  g(X)  in  the  case  of  normally  distributed  x2.  In 
this  case  the  lower  percentage  points  of  y^  are  (  ~  )  times  those  of 


When  £  is  not  small,  the  Pr(y1  ^  Q2),  (Q  ^  0),  must  be  considered. 
The  integrals  over  the  region  in  Figure  7  bounded  by  y^  =  Q2,  z  =  £  , 
yx  -  (1  +  <&(.)/G(z)  and  y^  =  ^2,  z  =  £  and  y^  =  (1  -  o0/G(z)  may  be 

12) 

A  detailed  account  of  these  integrals  is  also  given  by  Fisher 
(B.  A.  Tables,  1931). 
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expressed  as 


Pr(y1  4  q2) 


2o( 


-l/f' 


where: 


(f  -Y)  g(£  )df  - 


(f-X)g(f  )df 


Y  -  (1  +  ot  -  «i2)/(  >^q2),  X  =  (1  -  -  Q2)/(^Q21. 

4. 5. 2. 2.  The  Quotient  of  a  Triangular  and  a  Normal  Variable:  As  shown 
in  Section  4.2,  the  sum  (or  difference)  of  two  independent  rectangular 
variables  distributed  on  the  same  interval  is  triangularly  distributed. 
Considerable  difficulty  is  encountered  in  deriving  the  respective 
p.d.f. 's  of  Y.^  =  T/lN.  This  difficulty  is  attributable  to  the  fact  that 
the  p.d.f.  of  T  is  defined  in  two  branches  as  Equation  (4-39).  The  re¬ 
sulting  p.d.f.'s  are  functions  of  the  error  function.  In  the  simplest 
case,  i.e.JT  is  triangularly  distributed  on  (0,  2]  and  N  is  taken  to  be 
N(0,  l),the  p.d.f.  of  Y^  is  expressible  as  the  difference  of  two 


tabulated  error  functions.  So  that 

l/yl  2 

1  *y2 


[  ■ 

2/y, 


dy_  - 


0 

/ 

i/y. 


-fc  y- 


dy, 


71  /  0. 


( 4-50) 


Using  Broadbent's  method  of  analysis,  an  expression  convenient 
for  numerical  computation  may  be  derived. 

Consider  the  sum  (or  difference)  of  two  independently  distributed 
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rectangular  random  variables  defined  on  [a,  b]  with  the  same  half- 
range^^  •  Allow  the  sun  (or  difference)  of  their  means  to  be  de¬ 
noted  as  V  .  The  ratio  T/N  may  be  written  in  a  similar  form  to  that 
of  (4-43)  or 

y{  =  (  +  />  Y  )/( ^  +  <r  f  )  . 

(Conditions  being  suoh  that  K  >  0,  ^  =  e*  where  0<  2<*L  1,^  >  0f 

~  where  0  <  jA' <  oo).  The  variable  is  distributed  on 

[-2,  2]  and  }f  is  independent  of  Y£  and  satisfies  those  conditions 
outlined  in  the  preceding  section. 


The  standardized  quotient  is 

yl  (  X  yl  ~  x2 

II 

In  this  case, 

X  1  has  the  p.d. f. 

j  x1  -  1  +  2  *L 

;;.l 

^xx) 

4<*2 

/  1  ♦  2o<_  - 

4cxl2 

for  (1-2  e(<x1  ■£.  1)  t 


for  (1  -c  xx  -£1  +  2cxl  )  , 


The  joint  p.d.f.  of  z  =  F(*2)  and  y  *  is 

G(Z)  {  yx  G(Z)  -  1  +  2<  ]  dz  dy.jAcx',2,  ((.1  -  2c t  )/G(z))  4  y^l/GCz)) 

G(Z)  j  1  +  2oC  -  yx  G(Z)  |  dz  dy^io^2,  ((1/G(  z)  4  y1  ^  (1  +  2cxl  )/G(z)), 

when,  z  >€  and  corresponding  expressions  when  z  c  £  over  the 
region  shown  in  Figure  8. 


(4-5D 
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The  Pr(0  <  7l^  Qx).  >  O^is  the  integral  of  the  first 

element  over  the  region  A  bounded  by  ^  3  Qx*  z  -  l  and  y  =  (l-2cxl)/G(z) 
minus  the  integral  of  the  difference  of  the  two  elements  bounded  by 
3  Qp  z  -■  1  and  y^  =  1/G(z),  minus  the  integral  of  the  second 
element  over  the  region  bounded  by  y1  =  C^,  z  «=  1  and  y1  =  (l+2oO/Q(z) . 

The  first  of  these  integrals  is 


^2  f  *  l*.  1  [G(z)]  “  2Qj_  Q(»)  (1-2*0  +  (l-2o<J2jdz.  (4-52) 

z  =  F  jci-2oO/QjJ 

The  integral  is  simplified  by  changing  the  variable  of  integration 
to  £  and  writing 

U  -  (1  -  2<-  Q1)/(^Q1). 

The  integral  (4-52)  becomes 

v'V 

-T^r- 


Vz  j  (f  -  nr  g(  f )  d(  f). 

f=* 


Thus,  if 


00 

J2  g(D)  ij.  j~  (f  -  u)2  s(f>  d(f  ), 

f=° 

V  =  <l-Q1)/(  ^qx),and  V  -  (1+2«L  -Q1)/(  J^),  the  full  expression 
for  Pr(y1  c  Q  )  is 


8? 


It 


t£<2  2 

Pr(y.  6  q)  =  €  +  — — % 
1  (2o6)2 


J2  g(ll)  -  2  J2  g(V)  +  J2g(W) 


.  (4-53) 


may  be  concluded  hlso  that 


2 

Pr(y  >  Q)  =  - 

1  (2c^)2 


r  U2  +  1 


2^  "  J2  g(U)J  “  2  ^-~T^  ~  J2  b(V)  + 


Jo  s(W) 


-C. 


(4-54) 


Equations  (4-53)  and  (4-54)  have  similar  computational  advantages  as 
noted  for  (4-48)  and  (4-49) • 


V.  LIMITING  DISTRIBUTIONS,  APPROXIMATIONS  AND  ASYMPTOTIC  RESULTS 


5.1  Introduction 

A  frequently  encountered  problem  in  the  missile  and  space  industry 
may  be  described  as  follows: 

A  complex  mechanism  (say,  a.  missile)  contains  many 
individual  and,  to  some  extent,  independently  operating 
subsystems.  Such  structures  include  airborne  subsystems 
such  as  payload,  guidance,  electrical,  hydraulic,  flight 
control,  pressurization  and  airframe  as  well  aa  ground 
support  subsystems,  i.e.,  a  ground  power,  facilities  and 
propellent  subsystem. 

Frequently,  component  dataare  available  pertaining  to  failure  rates 
of  each  of  these  subsystems.  These  data  are  used  to  estimate  the  overall 
system  failure  rate.  Of  interest  then  is  the  statistical  statement 
which  may  be  made  concerning  overall  system  failure  probability. 

This  classical  problem  in  reliability  analysis  has  received  con¬ 
siderable  attention;  however,  most  attempts,  in  practice,  end  in  one 
of  two  methods  of  analysis.  There  may  be  a  decision  to  use  the  ex¬ 
ponential  failure  probability  density  function  to  estimate  the  failure 
probability  of  each  subsystem.  As  an  alternative,  the  use  of  a  digital 
computer  ia  made  to  simulate  operation  of  the  overall  system  and  thus 
gain,  through  repetition,  a  confidence  statement  of  the  overall  system 
failure  probability.  The  use  of  the  exponential  density  function 


features  the  highly  important  property  of  forming  a  complex  product 


simply  by  summing  the  exponents  of  the  variables  representing  com¬ 
ponent  failure  rates. 

In  general,  problems  concerning  the  establishment  of  a  reliability 
estimate  for  various  equipment  or  events  may  be  regarded  as  an  investi¬ 
gation  of  the  properties  of  a  random  variable  formed  from  the  product  of 
several  components.  Problems  of  this  type  may  be  investigated  by 
studying  the  distribution  of  Y^  which  may  be  expressed  in  the  general 
form 

h  -  n  hj 

As  seen  in  Chapter  IV,  Geary  [38]  produced  an  approximation  for 
Y1  =  N]/N2  wliere  N  represents  the  normal  density  function.  Fieller  [ 32 J 
shows  exact  percentages  of  N,/N_.  Craig  [  20]  and  Aroian  [2j  have  studied 
and  tabulated  results  of  In  addition,  Broadbent  [l4]  ,  [15]  studied 

measures  of  efficiency  which  are  of  the  form  R/N  and  T/N  where  R  and 
T  represent  the  Rectangular  and  Triangular  distributions.  Certain 
linear  combinations  of  random  variables  in  quotient  form  have  been 
examined  such  as  Gurland's  [48]  investigation  of 

K  =  (*l  X1  +  a2  X2  +  ,*,+  an  V/(bl  X1  +  b2  X2  +***+  bn  V 

where  the  joint  density  function  f(x^,  xp,  ,  x  )  is  known.  All 

of  these  results  are  of  a  specific  nature  and  are  not  generally  applicable 


n 

j=n+l 


(5-1) 


to  investigating  distributions  of  the  type  (5-1) • 

Several  approximations  to  distributions  of  random  variables  of  this 
type  (5-1)  have  been  investigated.  One  of  these  will  bow  be  described. 

5.2  Lognormal  Approximations. 

Consider  a  Measure  of  the  form  (5-1)  in  which  each  component  is 
subject  to  rectangular  or  normally  distributed  measurement  error  with 
known  half-range  or  coefficient  of  variation.  In  many  applications, 
the  distributions  of  the  sampling  or  measurement  errors  are  known  to  be 
one  of  these  two  forms,  and  their  standard  deviations  or  half-ranges 
can  usually  be  found  by  simple  investigation.  In  other  applications, 
the  distribution  of  these  errors  can  only  be  estimated  and  then  per¬ 
haps  the  roughest  approximations  should  be  used  only  as  bounds  for 
experimental  error. 

It  is  impossible  to  list  exhaustively  all  possible  combinations 
of  errors  to  consider  in  (5-1) •  In  the  event  of  a  large  number  of 
errors  being  combined,  the  asymptotic  distribution  may  be  used  with 
confidence  Jcramer  [221]].  It  is  in  the  intermediate  cases  that  ap¬ 
proximations  must  be  crii.ically  considered.  The  choice  of  a  suitable 
family  of  approximating  distributions , however,  will  always  be 
rather  subjective.  In  this  respect,  the  consideration  of  an 
approximation  is  directed  towards  finding  a  "working,  first-order 
approximation".  Several  authors,  i.e.,  Shellard  [115 J  »  Finney  [55), 
Gaddum,and  Johnson,  have  investigated  the  use  of  lognormal  approximations 
to  (5-1)  as  an  alternative  to  attempting  to  find  its  exact  distribution 
using  the  methods  of  Section  3.^.  They  have  shown  that  the  p.d.f.  of 


Y^  (5-1)  tends  to  the  lognormal  as  s  - ►  oo  under  very  general 

conditions. 

The  most  widely  used  lognormal  approximation  to  (5-!)  is  *  variate 

z  such  that  log(z-  f*  )  is  normally  distributed  with  mean  /i  and  variance 

2  2 
<J~  .  Here,  for  simplicity,  only  the  methods  of  selecting  and  cr 

are  considered,  f  being  considered  as  equal  to  zero.  The  notation 


oC^,  •  »  •  ,  odn  denotes  the  coefficient  of  variation  or  the 

quotient  of  the  half -range  and  the  mean  of  a  rectangular  variable  for 


each  of  the  variates  x^,  x^, 


•  *  i 


x  .  Thus  Y..  =  R/N, 

n 


(<  =°  ♦  33,«12  =  0.10),  denotes  the  quotient  of  a  rectangular  variate, 
whose  mean  is  three  times  its  half-range  and  an  independent  normal 
variate  whose  coefficient  of  variation  is  10#.  The  variate  z  is 
necessarily  positive  although  may  have  negative  values.  When 
o^l ,  o^,  .11,  are  small  positive  values,  the  approximation  by  this 
lognormal  distribution  may  be  satisfactory,  since  the  probability  attached 
to  the  possibility  of  negative  values  is  exceedingly  small. 

There  are  two  methods  of  fitting  the  lognormal  approximation.  The 

2 

first  method  of  choosing yi  and  cT*  is  to  calculate  the  moments  of 

2 

log  y^  and  to  set yi  and  <T  equal  to  the  first  and  second  corrected 
moments.  This  method,  termed  the  method  of  fitting  by  moments  to  log  y^, 
was  first  detailed  by  Finney  [34]  (19^1)*  An  alternate  method  is  to 
choose  the  lognormal  approximation  whose  mean  and  variance  are  equal 
to  the  mean  and  variance  of  y^.  This  method  generally  termed,  the 
method  of  fitting  by  moments  to  y  ,  differs  from  that  of  the  first 


method  but  the  difference  is  seen  to  disappear  as  n  increases#  The 
possibility  of  two  methods  raises  the  question  of  which  method  should  be 
used. 

2 

Let  a  lognormal  p.d.f.  have  mean  m  and  variance  v  and  lot  the  p.d.f. 

of  the  logarithm  of  the  variate  from  the  lognormal  have  mean  ;u  and 
2 

variance  cr  .  it  is  easily  shown  that 

2 

m  =  exp  ( p  +  a~  /2) , 

v2  =  exp  | (2^  +  <r2)}  fexp(  or*2)  -  l}  ,  (5-2) 

2  2 

p  =  log  m  -  #  log  (1  +  v  /«  ), 
cr2  =  log  (l  +  v2/m2). 

The  fitting  to  log  y1  requires  that  we  select  the  mean  and  variance  of 

2 

log  y  and  set  p  and  <r  equal  to  them.  The  k#  point  of  this  fit  is 
exp( p  +  p^  cr),  where  the  probability  that  a  standardized  normal  variable 
is  less  than  is  t/100. 

The  fit  to  y1  requires  that  the  mean  and  variance  of  y^  be  found 
2 

and  set  equal  to  m  and  v  .  The  k#  point  of  this  fit  is 

exp  [pt  |  log  (1  +  t^)}  ]J/  (1  +  t2)  ,  (5-3) 

where:  t  =  s/m,  the  coefficient  of  variation  of  y^. 

The  k%  points  of  (5-3)  are  tabulated  in  Table  VIII,  Appendix  A.  The 
appropriate  points  are  found  by:  (i)  determining  the  expected  value  and 
coefficient  of  variation  of  y^,  (2)  entering  the  table  with  the  appropriate 
t,  and  (3)  multiplying  the  value  in  the  table  by  m. 


In  order  to  find  the  momenta  of  or  of  log  y^  expressions  are  derived  for 
the  moments  of  various  positive  and  negative  powers  of  the  normally 
and  rectangularly  distributed  variables.  These  variables  may  be  ex¬ 
pressed  as  x  =  ji(l  +  o/w)  where  o(.and  w  are  as  previously  defined. 

In  the  case  of  normally  distributed  w,  the  moments  of  log(l  +  olw) 
are  more  easily  derived  from  the  truncated  distribution, 


exp  <-fc  w2)  d«/  (d-O/  JTn  j 


y  >  -  -L.  (5-4) 

oC  * 


which  has  the  property, 

Pr(l  +  o<w  <■  0)  =  €  ,  (€  <  lCf^  when  0  <  o(i  0.15)* 

For  all  practical  purposes,  (5-^)  ie  indistinguishable  from  the 
normal  p.d.f.  with  the  above  restrictions. 

The  raoraent-geherating  function  of  log(l  +  ®^w)  ie  E  ((1  +  o(.w)  J 
and  may  be  written  as 


*  r 

j2T(l-€)  J 


(1+it  o(.w  +  Yz  it(it-l)o<L2  w2+. .  .)exp(-#  w2)dw  +  K  (5-5) 


where  |K|^  €  /(l-€). 

Broadbent  [l4]  has  discussed  this  integral  and  its  convergence 
and  has  ahown  that  it  leads  to  the  cumulants  of  log  (1  +  o/w).  He 
has  also  shown  that  the  cumulants  of  log  (1  +  c/w)  when  w  is  rec¬ 
tangularly  distributed  may  be  obtained  and  that  a  similar  analysis  leads 


wmm 


mmm 


to  the  expected  values  of  various  powers  of  (1  +  o£w)  for  both  cases 
13 

of  w.  These  expressions  for  E  |(1  +  <^w)rj  and  E  [log(l  +  <Xw)r] 

for  r  =  +)4,  +  1,  +  2,  and  +  4  are  shown  in  Table  IX,  Appendix  A. 

^®^ore  considering  a  complicated  case  for  the  lognormal  approxi- 
mation,  it  is  worthwhile  to  consider  the  very  simple  independent  case 
Nx  (1  +  oL1  w) 

of  Y1  =  Nj  =  "(1"+'  oc2  «)  where represent  the  normal  p.d.f.  The 

coefficidmts  of  variation,  0^  and  o^2  of  the  numerator  and  denominator 
are  taken  as  100  and  100  c*^  respectively. 

Using  Table  IX,  these  results  are  obtained: 

2  H = E  H  E  ft] '  ■ =  (i>  a  *  °l*  * 3  +  •••>. 

E  [r/j  =  E  [y]  *-3]  =  <u  ^)(W  15  ^ 

V  -  *|V]  -  E  2  -,2  .  ax2  +  oc/oo^2  «j2*8  . 

Similarly  for  the  method  of  fitting  to  log  y  ,  we  obtain  these 


results: 


E  |"lo6  yxj  =/»  =  -Cfc  ^  0^^+...)  +  (#  ^  +  ...)  , 

V  [log  yj  =  v2  =  (od,2  +  |  06^  +...)  +  (o^2  +  |  od24  +  ...). 


J.  Haldane  j*5i  ]  has  also  derived  these  moments. 


The  first  two  moments  of  the  lognormal  fitted  to  log  y^,  using  these 
2 

values  of  )i  and  cr  ,  are: 

Mean  =  1  +  o oL ^  |  oL^  +  ...  , 

Variance  =  o^2  +  +  3  +  3  e^2  o(.^  +  5  +  . .. 

2  2  2 

These  estimates  agree  with  m  and  v  to  0(o£,^  ,  ot2  ),  and  to  higher 
ordera  when  o<^  =  o^2»  thus  the  differense  in  the  two  fits  is  small  when 
oC^  and  oC^  are  small.  In  fact,  for  this  case,  there  is  complete  agree- 
ment  in  the  first  three  significant  digits  at  the  1$6,  3%  and  the  95#  points 
between  the  two  methods  of  fitting  and  the  exact  points  computed  by 
Craig's  methods  (4-8).  The  difference  between  exact  points  and  the 
two  methods  of  lognormal  'Approximation  may  be  shown  to  be  ^  +  ,002 
at  the  1$,  5#»  95#  and  99#  points  in  the  cases  of  R^/EL,i  RN,  and  N^/N,,. 

We  now  return  to  the  problem  of  setting  limits  to 

'  ui  *2  •  •  •  xi)Axyi  ■  •  •  *»>• 

If  the  number  of  component  variates  is  larger  than  two,  the  log¬ 
normal  approximation  vdll  give  satisfactory  results  when  the  coefficients 
of  variation  of  the  components  are  small  and  do  not  vary  greatly.  It 
is  , therefore,  a  suitable  approximation  for  general  use. 

To  calculate  the  first  two  moments  of  y^,  it  is  necessary  to  know 
the  coefficient  of  variation  or  quotient  of  half-range  and  mean,  of  each 


VI.  CHARACTERIZING  PROPERTIES  OF  STATISTICAL  DISTRIBUTIONS 


6.1  Introduction. 

If  X^  and  X^  are  independent  and  +  X„  is  known  to  obey 

a  normal  p.d.f.,  it  is  well  known  (Cramer  [22]  )  that  and  must 
each  be  normally  distributed.  This  suggests  the  general  problem  of 
determining  when  the  p.d.f.  of  an  algebraic  combination  of  random 
variables  uniquely  determines  the  respective  p.d.f.' a  of  the  component 
random  variables.  Investigations  of  this  type  have  led  to  a  class  of 
statistical  problems  formally  termed  "characterizing  properties  of 
statistical  distributions." 

An  important  phase  of  this  problem  is  that  of  determining  the 
properties  of  observations  from  their  estimated  p.d.f.  when  it  is  known 
that  the  samples  are  in  reality  either  a  product  or  a  quotient  of 
random  variables. 


6.2  Three  Important  Problems. 


By  Theorem  3t  Chapterlll,  the  p.d.f.  of  a  random  variable  quotient 
Y^  -  X^/X_, ,  where  X^  and  are  arbitrary  independent  random  variables, 
is  given  in  its  most  general  form  as 


where  the  transformation  y^  =  x^/x^,  y^  ~  x~>  implied  and  (  *  ) 
and  ( •  )  are  the  respective  independent  p.d.f's  of  X^  and  Xp. 


Several  important  problems  may  be  formulated  regarding  (6-1); 

three  of  which  will  now  be  briefly  considered . 

The  first  and  most  straight-forward  is  that  of  determining  ^(y^)  when 

the  functional  forms  of  (  ( *  )  and  gY  (  •  )  are  known.  This  prob- 

X1  2 

lem  simply  reduces  to  the  integration  of  (6-1). 

The  second  problem  which  many  consider  as  being  the  "classical” 
characterization  problem  may  be  formulated  as  follows: 

Let  X^  and  X.,  be  two  independent  and  identically  distributed 
random  variables  having  a  common  p.d.f.,  f(x).  Let  the  quotient 
^1  *  ^S/^2  be  ^ stributed  insome  known  p.d.f.,  ^(y^).  Is  it  possible 
to  obtain  a  characterization  of  f(x)  by  the  properties  of  the  quotient 
p.d.f? 

In  this  case,  Equation  (6-1)  may  be  written  as 

/oo 

f(yl  y2^  f(y25  I  ^2 1  d>2  (6-2) 

-oo 

where  only  the  form  of  y^ )  is  known.  Among  the  questions  which 

naturally  arise  are:  Cl)  how  many  independent  solutions  of  f(*  ) 
satisfy  the  linear  integral  equation  (6-2)  for  a  given  function  of 
'^1  )^’  and  if  ,nore  than  one,  is  it  possible  to  deduce  the 
general  properties  of  the  class  of  c\H  f(x)  satisfying  (6-2)? 

Several  investigators  have  dealt  with  a  Cauchy  form  for  $?(y^). 

When  X^  and  are  independently  and  normally  distributed  random 
variables  with  parameters  (0,  1)  or  ^0,  a*  ),  it  is  well  known  that 


is  Cauchy  distributed  symmetrically  about  y^  =  0.  It  has 

been  conjectured  that  this  is  a  unique  property  of  the  normal  p.d.f.; 
a  hypothesis  which  implies  the  existence  of  only  one  solution  of  f(x) , 
namely,  f(x)  =  N(0,  <T  ),  which  satisfies  the  linear  integral  equation, 


x  tr  y 


OO 


‘(y2)f(yIy2)  |  y2 1  dy2»  (C  =  a  constant),  (6-3) 


-oo 


Laha  [75],f(6],  StecViJli7j  and  Mauldon  j8^  all  have  subsequently 
proved  this  conjecture  to  be.  false  by  exhibiting  several  non-normal 
solutions.  Goodspeed  [4oj  has  shown  in  fact  that  an  infinity  of 
solutions  to  (6-3)  exists. 

Laha  (?6j  has  shown  that  if  and  X2  are  independently  and 
identically  distributed  random  variables  having  a  common  o.d.f. 

F^x)  such  that  the  quotient  Y^=X^/X2  follows  the  Cauchy  Law  distrib¬ 
uted  symmetrically  about  zero ,  then  F(x)  has  these  gsncral  properties! 

l)  F(x)  is  symmetric  about  the  origin  x  =  0  ; 
ii)  F(x)  is  absolutely  continuous  and  has  a  continuous  p.d.f. 
f(x)  s  F'  (x)  >  O  5 

iii)  the  random  variable  X  has  an  unbounded  range; 
iv)  the  p.d.f.  f(x)  satisfies  the  linear  integral  equation 


oo 


1+  y. 


f(yl  y2^  f^y2l  y2  dy2  ’ 


(6-4) 


Mb 


100 


mmw 


PUSMSI 


H  nrU  «t  U  .  -  rin.,1  I  j 


1111 


which  holds  for  all  y^,  where  C'  is  a  constant.  (Here  again  the 
transformation,  y^  =  x^/x^,  y2  =  x 2  is  implied). 


It  is  now  apparent  that  an  answer  to  the  question  of  the  number  of 
f(x)'s  satisfying  the  Qquation  (6-4)  is  the  equivalent  of  completely 
enumerating  all  solutions  to  this  integral  equation  This  problem  is 
very  difficult  and  still  remains  to  be  solved. 

Mauldou(T35  3  has  studied  similar  problems  showing  in  particular 
that:(l)  the  ratio  of  two  independently  and  identically  distributed 
chi-square  variables  is  not  the  only  quotient  distributed  according  to 
the  general  F-distribution, 


^(yx) 


C  (a+b) 

r(a)Rb) 


b-1 

77— ».  b>°) 


(6-5) 


and  (2)  that  we  cannot  assert  that  f(x)  is  normal  even  when  a  sample 


of  x^,  i  =  1,  .  .  .  ,  n  observations  is  taken  from  the  same  p.d.f. 
distributed  symmetrically  about  zero  such  that  the  cample  statistic 
t-l  yvT  is  distributed  as  Student’s  t-distribution  (x,  a 
denoting  the  usual  sampling  statistics.) 

The  third  problem  is  a  general  case  of  the  preceding  problem  and 
involves  solving  the  integral  equation 


00 

w  -f  vyi 


y2)  gx2(y2)  |y2l  *y~ 


(6'ij 


-00 


when  one  of  the  two  combinations  of  functions,  (f^Xy^),  f ^  (  •  ))  or 
(S^?(y1) ,  gx  (  •  )),  is  known. 
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In  certain  sampling  studies,  it  is  conceivable  that  collected 

data  y. ,(i  =  1,  .  .  .  ,  n),will  be  of  a  ratio  form,  say  y.  =  x  /x  , 

1  1  di 

and  at  the  same  time  estimates  of  the  respective  p.d.f.'s  of  the  ratio 
y^  and  of  one  of  the  components,  say  can  be  made  with  a  high 
degree  of  certainty.  In  these  cases,  it  is  reasonable  to  question 
whether  or  not  the  general  properties  of  X^  may  be  ascertained. 

A  simple  method  of  obtaining  an  estimate  of  the  p.d.f.  of  X.,  is 
to  consider  the  ratio  =  X^/Y^  using  the  estimated  p.d.f.  of  the 
respective  components.  However,  the  p.d.  f .  (  )  may  reprebent 

only  one  of  several  possible  solutions.  In  this  case  only  from  the 
complete  enumeration  of  all  solutions  satisfying  equation  (6-6)  will  it 
be  possible  to  deduce  the  general  properties  of  gy  (  •  ) . 
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TABLE  VI il 

Standardized  Lognormal  Percentage  Points 

Given  the  mean  o  and  standard  deviation  ,  let  t  =  100  s/a.  The  per¬ 
centage  point  of  the  lognormal  distribution  with  this  mean  and  standard 
deviation  is  the  entry  in  the  table,  multiplied  by  a* 


0  1.0000 

1  0.9770 

2  0.9544 

3  0.9322 

4  0.9105 

5  O.M892 

6  0.8633 

7  0.8478 

8  0.8278 
9  0.8082 

10  0.//890 

11  0.7702 

12  0.7518 

13  0.7338 

14  0.7162 

15  C.6900 


1.0000 

0.9836 

0.9675 

0.9514 

0.9356 

0.9200 

0.9045 

O.K892 

0.8741 

0.9592 

0.8445 

0.8299 

0.8156 

0.8015 

0.7875 

0.7738 


Lower  1#  points  (l 

C  =  1) 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

9954 

9930 

9907 

9884 

9861 

9838 

9724 

9701 

9679 

9656 

9633 

9611 

9499 

9477 

9454 

9432 

9410 

9388 

9278 

9256 

9233 

9213 

9169 

9148 

9061 

9040 

9019 

8997 

8976 

8955 

8849 

8828 

8807 

87 8? 

8766 

8745 

8641 

8621 

8600 

8580 

8559 

8539 

8438 

8418 

8398 

8378 

8357 

8338 

8238 

8219 

■8199 

8179 

8160 

8l4o 

8043 

8024 

8oo4 

5985 

7966 

7947 

7852 

7833 

7814 

7795 

7776 

7758 

7665 

7646 

7628 

7609 

7591 

75 72 

7481 

7462 

7445 

742? 

7409 

7391 

7302 

7285 

7267 

7249 

7232 

7214 

7127 

7110 

7093 

7075 

7058 

7041 

Lower  5%  points  (k  =  5) 


0.3  0.4  0.5  0.6 


9983  9967 
9820  9804 
9658  9642 
9498  9482 
9340  0324 
9184  9168 

9029  9014 
8877  8862 
8726  8711 
8577  8562 
8430  8415 

8235  02'?1 
8142  8128 

8001  798? 
7861  7848 


0.7  0.8 


9885 

9723 

9562 

9403 

9246 

9091 

8937 

8786 

8636 

8489 

8343 


9853 

9691 

9530 

9372 

9215 

9060 

8907 

8756 

8607 

8460 

8314 


8242 

8099 
7959 
78201  7806 


t  o.o 


Upper  1%  pointed  =  99) 


0.2  0,3  0*4  0.5  I  >.6  0.7  0.8 


1.0000 


1606 
1870 
2111  2138 
2384  2411 
2661  2689 

2942  2970 
3228  3257 
3518  3547 
3813  3843 


1659 
1923 
2165  2193 
2439  2466  2494 
2716  2?45  2773 

2999  3027  3056 
3285  3314  3343 
3577  3606  3636 
3872  3902  3932 


Upper  5* 


0.1  0.2  0.3  0.4  0,5  0.6  0.7  0.8  0.9  A 


1.00001  0016  0033  0049  0066 

0182  0198  0215  0232 

0349  0366  0383  0400  041? 

0518  0535  0552  0569  O586 

0688  0705  0723  0740  0757 

0860  0878  0895  0912  0930 

1034  1051  1069  1086  1104 

1209  1226  1244  126 2  1279 

1385  1403  1421  1439  1456 

1563  1581  1599  1617  1635 

1742  1760  1779  1797  1815 

1923  1941  1959  1978  1996 

2104  2123  2142  2160  2178 

2288  2307  2325  2343  2362 

24?2  2491  2510  2528  2547 

15  I  1,26391  - 


0115  0132 
0282  0299 
0450  0467 
0603  0620  0637 
0774  0791  O809 
0947  0964  0982 

1121  1139  1156  |  1174 
1297  1315  1332 
1474  1492  1510 
1653  1671  I689 
1833  1851  1869 

2014  2C32  2050 

2196  2215  2233 
2380  2399  24l? 

2565  2584  2602  |  2621  |  +19 


TABLE  X 


Tables  of  the  cumulative  distribution  funot  on  of  the  random 


variable  Z  -  — -  ,  where  X,  are  normally  distributed. 

*i  % 

The  p.d.f,  of  Z  ie  f<a)  =  IjU)  -  ^(z)  ia  *hich 


♦ 

and  I2(z)  is  the  same  integral  defined  on  the  interval  (-00,  0). 


Values  of  f(z),  the  probability  density  function,  andF(z),  the  cumulative 


distribution  function  of  z  for  various  values  of  v. 


and  the  coefficient  of  correlation  p  . 


TABLE  X.a 


Case  1:  ^  a  0,  =  0,  vp  s  0.  The  p.d.f.  of  f(z)  ia  symmetrically 

distributed  about  the  singularity  point  z  -  0  with  U(z)  =0,  c r  =  1. 

z 

In  tide  case  F(z^  =  1  -  F(-z). 


.00005 
.00011 
.00026 
.00061 
.00146 
.00355 
.00557 
.00878 
..01.395 
.02235 
. 03625 
.04645 
.05983 
.07756 
.10138 
.13402 
.15493 
.17996 
.21025 
.24749 
.29425 
.35477 
.43887 
.55790 
.77256 
+0* 


.00004 
.00019 
.00021 
.00056 
.00139 
.00323 
.00509 
.00796 
.01231 
.01943 
.03090 
.03914 
.04969 
.06336 
.08114 
.10450 
.11890 
.13144 
.15508 
.17790 
.20491 
.23721 
.27 656 
.32585 
,45242 
.50000 


Mi 


TABLE  Xb 


'ffl 


Case  2:  ^  -  0,  ^  *  (1.0)(0),  v2  =  0(1.0).  The  p.d.f.,  f(z)  is 
symmetrically  distributed  about  the  singularity  point  z  =  0  with 
E(z)  =  0  and  <x %  =  "\fz>  in  oao*  = 


2 

f  (z) 

F(  z) 

-12.0 

.00001 

.00001 

-11.0 

,00002 

.00002 

-10.0 

.00005 

.00004 

-  9.0 

.00012 

.00013 

-  8.0 

.00029 

.00034 

-  7.2 

.00058 

.00069 

-  6.4 

.00118 

.00139 

-  5.6 

.00239 

.00282 

-  4.8 

.00485 

.00576 

-  4.0 

.00981 

.01166 

-  3.6 

.01397 

.01642 

-  3.2 

.01992 

.02319 

-  2.8 

.02846 

.03287 

-  2.4 

.04078 

.04572 

-  2.0 

.05873 

.06662 

-  1.8 

.07069 

.07949 

-  1.6 

.08533 

.09509 

-  1.4 

.10340 

.11397 

-  1.2 

.12595 

.13690 

-  1.0 

.15460 

.16490 

-  0.9 

.17195 

.18129 

-  0.8 

.19193 

,19948 

-  0.7 

.21519 
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The  components  of  a  quotient  or  product  of  a  random  variable 
assume  many  different  representations.  Not  only  ia  one  interested 
in  ratios  and  products  of  random  variables  described  by  a  specific 
frequency  function,  but  also  of  these  functional  forms  of  variances, 
ranges,  proportions,  various  forms  of  Milla*  ratio,  correlation  co¬ 
efficients,  etc. 

The  following  bibliography  is  believed  to  be  very  comprehensive 
in  respect  to  articles  pertaining  to  the  distribution  theory  of 
products  and  quotients  of  random  variables.  However  only  a  small 
sample  of  articles  pertaining  to  other  subjects  as  those  mentioned 
above  are  given. 

An  index  is  presented  in  which  listings  appear  by  headings 
classifying  many  of  the  various  combinational  forms  of  quotients 


and  products 
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1.  Aroian,  L.  A.,  "A  Study  of  R.  A.  Fisher's  Z  Distribution  and  the 
Related  F  Distribution,"  Annals  of  Math.  Stat.  Vol.  12,  19^1. 

The  purpose  of  this  paper  is  to  di souse  the  semi-invariants  of  the  Z 


distribution  and  to  fin^  useful  approximations  for  them. 

B1  2 

bution  of  Z  =  Yz  log  — g  , where  ,(i  =  1,  2), are  sample 
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variances  ,1s  well  known; 


p(z) 


tfn,  yn- 

1  £ 

2nl  n2 


n^z 


e2z  +  \ 


The  distri- 
estimates  of 


dz. 


The  author  shows  that  as  and  n^  approach  infinity  in  any  manner  the 
distribution  of  Z  approaches  normality. 


2.  Aroian,  L.  A.,  "The  Probability  Function  of  a  Product  of  Two 

Normally  Distribution  Variables,"  Annals  of  Math,  Stat.  Vol.  18,  19^7. 

Let  X^  and  X^  follow  a  normal  bivariate  probability  function  with 

means  Xn  and  X_  standard  deviations  <F.  and  01  and  correlation  co- 

efficient  r  ^  ,  Also,  let  ^  and  =  X^/o^.  Craig  jsij 

found  the  probability  function  of  Z  =  X^  X^/  0^  o^  in  closed  form  as 
the  difference  of  two  integrals.  Craig,  for  the  purposes  of  numerical 
computation,  expanded  this  result  in  an  infinite  series  involving  powers 
of  Z,  and  Bessel  functions  of  a  certain  type.  Difficulty  arises, 

for  large  and  p,,,  the  convergence  of  this  expansion.  Aroian 
demonstrates  that  Z  may  bo  approximated  by  the  Gram-Charller  series 
and  the  Type  III  function  and  presents  the  percentage  points  of  F(s)  for 
special  case  =  1,  V^  =  10  and  r  =  0. 
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3.  Aroian,  L.A.,  Tables  and  Percentage  Points  of  the  Distribution. 
Function  of  a  Product,  Hughes  Aircraft  Company,  California,  1957. 

Let  X^  and  X 2  be  normal  uncorrelated  variables  with  means  and 

m2  and  variances  and  tr,,  respectively,  and  let  =  m^/o^ .  The 

author  lmdsi  (1)  The  100  oC percentage  points  of  the  distribution  of 

2  =  XiX2//,°i°2  are  tat)Uj-ated  for  different  values  of  S  (i  =  1,2);  (2) 

The  cumulative  distribution  of  Z  is  also  tabulated. 

.  it 

4.  Baker,  G.  A.,  Distribution  of  the  Means  Divided  by  t^e  Standard 

Deviations  in  Samples  From  Nonhomogeneous  Populations  Annals  of 
Math.  Stat.  Vol.  3,  .1.932.  - 

It  is  the  purpose  of  this  paper  to  discuss  the  distribution  of  the 

means  of  samples(of  size  two)  measured  as  the  mean  of  the  population 

which  have  been  divided  by  the  standard  deviations  of  the  samples. 

Experimental  results  for  samples  of  four  from  normal  populations  are 

presented, 

tt 

5.  Baker,  G.  A.,  Distribution  of  the  Ratio  of  Sample  Range  to 
Sample  Standard  Deviation  for  Normal  and  Combinations  of 
Normal  Distributions,"  Annals  of  Math.  Stat.  Vol.  17,  ±946. 

This  note  summarizes,  in  tabular  form,  some  of  the  numerical  results 

obtained  m  previous  studies  of  the  distribution  of  sample  ranges  m 

terms  of  the  standard  deviation  of  the  sampled  population  for  homogeneous 

populations, 

fl 

b.  Banerjee,  S,  K. ,  A  Lower  Bound  to  the  Probability  of  Student's 
Ratio t  Sankhya ,  Vol.  18,  1957. 

Samples  of  size  n  are  drawn  from  a  population  having  mean  m,  and 
2nd  and  4th  cumulant3  k?  and  k^  respectively.  Starting  with  Tcheby cheff 's 
lemma  that  if  a  variate  X  has  mean  ra  and  variance  <r*  then 

P(X  6  m-k2)  <£  — — 
cr >  kH 
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whatever  W~  may  be,  the  behavior  of  Student's  ratio  i3  considered  from 
the  viewpoint  of  estimating  the  confidence  interval.  In  reticular,  a 

lower  bound  to  the  probability  o.f  the  event  x  +  s:  m la 

\jn  \n 

1  k4 

calculated  to  be  -  >  where  B0  -  — p  t-  3. 

V5  ,2  |  t'  ‘  k2 

»  1 n-1  ! 

A  table  is  given  for  t  =  3  and  several  values  of  and  n. 
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7.  Barrow,  D.  T.,  and  Cohen,  A.  C.  Jr.,  "On  Some  Functions  Involving 
Mills'  Ratio,"  Annals  of  Math.  Stat. .  Vol.  25,  1951*. 

In  this  note,  it  is  proved  that  for  all  finite  values  of  h, 
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te 
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*(h)  =  X . 

(Z-h  T 

2 

is  monotonic  increasing  so  that  2  >  0  and  1  «=  \j/{ h)  *e  2.  The 

function  Z  is  the  recipiocal  of  Mills'  ratio  and  and  m ^  are  the  first 
and  second  moments  of  a  singly  truncated  normal  population  about  a  point 
of  truncation.  The  function  'l'(h)  frequently  arises  in  connection  with 
maximum  likelihood  estimates  of  population  parameters  from  singly  trun¬ 
cated  normal  samples. 

8,  Bartlett,  H. ,  "The  Distribution  of  Second  Order  Moment  Statistics 
in  the  Normal  System,"  Proc.  Camb.  Phil.  Soc.,  Vol.  28,  1932* 

Let  x  be  a  normally  distributed  variable  and  the  distribution  is 

written  as 

d  p  =  e_!4x  /k2  dx 

where  k_,  is  the  serai -invariant  of  order  2.  The  moment  generating 

2 

function  of  the  distribution  of  x  is  developed,.  From  this  ,the  generating 

2 

function  of  the  semi-invariants  (k  )  of  x  is  developed. 
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9.  Barton,  D.  E. ,  David,  F.  N. ,  and  O'Neill,  A.,  Some  Properties  of  the 
Logarithm  of  the  Non-Centrai  F  ,  Biometrika,  Voi.  47,  J-960. 

2 

Given  two  independent  variables:  6  distributed  an  u  non-central 

chi  square  variable  with  f^  degrees  of  freedom  and  non-central  parameter 
2 

and  X  distributed  as  a  control  chx  square  variable  with  f'2  degrees 
of  freedom,  the  non-central  F  is  defined  as 

u2/f.) 

(x2/f2) 

The  log  transformation  of  F  and  the  possibility  of  a  normal  approx- 
amation  to  the  transformation  are  discussed  here. 

it 

10.  Basu,  D.  On  the  Dependence  of  Linear  Functions  of  Independent 
Chance  Variables,"  Bull.  Inter.  Stat.  Inst.,  Vol.  33,  1953. 

Relying  on  the  assumption  that  the  random  variables  X,,X_,  »...,X 

l  c.  n 

are  independently  distributed  and  possess  finite  moments  of  all  order, 

the  author  proves  the  following:  If  there  exist  two  linear  functions 
n  n 

f  =  Xj  a.  x,  and  t  =  X>  b.  X.  with  a  b  /  0  and  fund  £ 
j=l  j»l  5  J  3 

are  stochastically  independent,  then  each  x.  is  normally  distributed. 

3 

n  ii 

11.  Bliss,  C.  I.,  The  Calculation  of  Dosage-Mortality  Data  ,  Annals 
of  Applied  Biology,  Vol.  22,  1935» 

Very  detailed  applications  of  experimental  data  to  the  problem 
of  estimating  the  quotient  of  two  quantities  subject  to  experimental  error 
are  given. 

fl 

12.  Bose,  S. ,  On  the  Distribution  of  the  Ratio  of  Variances  of  Two  ,, 
Samples  Drawn  From  a  Given  Normal  Bivariate  Correlated  Population, 
Sankhya,  Vol.  2,  1935* 

A  formulu  for  the  p.d.f.  of  the  bivariate  chi  distribution  is  given 
This  formula  is  corrected  by  Knshnniah,  et.  ai.  [?o].  Tne 


distribution  of  Z  =  U/V,  the  ratio  of  two  correlated  chi  variates  io 
given. 

13.  Boyd,  A.  V. ,  "Bounds  for  Mills'  Rutio  for  the  Type  III  Distribution," 
Annals  of  Math.  Stat.,  Vol.  29,  1958. 

Des  Raj  [ lOl]  and  Cohen  have  shown  that  in  estimating  the  parameters 
of  the  truncated  Type  III  populations,  it  is  necessary  to  calculate  for 
several  values  of  x,  the  Mills'  ratio  of  the  ordinate  of  the  standardized 
Type  III  curve  at  x  to  the  area  under  the  curve  from  x  to  00 .  For  large 
values  of  x,  existing  tables  are  inadequate.  The  object  of  this  note  i/s 
to  establish  lower  and  upper  bounds  for  this  ratio. 

14.  Broadbent,  S.  R.,  "Lognormal  Approximations  to  Products  and  Quotients," 
Biometrika,  Vol.  43,  1956. 

The  problem  considered  is  to  obtain  bounds  of  limits  upon  combinations 
of  random  variables  (random  in  the  sense  of  measurement  error)  in  a  product 
or  quotient  form.  The  general  form  is  considered  to  be 

q  =  (x,  •  x„  .  .  .  x)/(x.«x  .  .  .  x  )  ,1  f  n  <  r  . 

12  n  n+1  m  r  '  — 

A  summary  of  the  known  exact  results  is  first  given.  The  lognormal 

distribution  which  is  asymptotically  exact  is  shown  to  give  useful 

approximations  when  fitted  by  moments  to  the  combination  above.  Tables 

arc  given  which  make  its  application  relatively  simple. 
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31.  Epstein,  Benjamin,  Some  Applications  of  the  Mellin  Transformation 
in  Statistics’*  Annals  of  Math.  Stat.  ,  Vol.  19,  1948„ 

It  is  well  known  that  the  Fourier  transform  is  a  useful  analytic 
tool  for  studying  the  distribution  of  the  sum  of  independent  random 
variables.  It  is  the  purpose  of  this  paper  to  study  the  Mellin  Trans¬ 
form  in  relation  to  products  of  independent  random  variables  <~.f  the  form 
V  =  i where  y  is  defined  on  positive  intervals.  In  this  paper 
Epstein  uses  the  Mellin  transform  to  obtain  V  =  N^N^  as  a  Bessel  function. 

it  n 

32.  Fleller,  E.  C.,  The  Biological  Standardization  of  Insulin, 

Journal  of  the  Royal  Stat.  Soc,,  (Series  B),  Vol.  7,  1940. 

Results  involving  quotients  of  random  variables  and  the  p.d.f.  of 

an  index  are  used  in  an  application  of  a  biological  assay  of  insulin. 

it 

33*  Fieller,  E.  0.,  The  Distribution  of  an  Index  of  a  Normal  Bivariate 

ii 

Population,  Biometrika,  Vol.  20,  1932. 

This  is  a  notable  paper  concerning  the  ratio  of  two  normally  distrib¬ 
uted  random  variables.  By  quadrature,  Fieller  presents  a  solution  to 
the  p.d.f.  of  =  N^/N^  for  the  general  case.  The  normal  case  is  expressible 
in  terms  of  the  bivariate  normal  p.d.f. 

34.  Finney,  D.J. ,  "The  Distribution  of  the  Ratio  of  Estimates  of  Two 
Variances  in  a  Sample  From  a  Normal  Eivariate  Population, 

Biometrika,  Vol.  30,  1938. 

Finney  shows  the  p.d.f.  of  the  ratio  Z  =U/V  of  two  correlated  chi  variates. 
From  this  result,  he  obtains  the  c.d.f.  of  Z  and  discusses  certain  pertinent 
applications. 

If 

35.  Finney,  D.  J.,  On  th^  Distribution  of  a  Variate  Whose  Logarithm  is 
Normally  Distributed,  Journal  of  the  Royal  Stat,  Soc.,  ( B) ,  Vol.  ?,  1940. 

The  problem  considered  here  is  to  determine  efficient  estimates  of 


both  the  mean  and  the  variance  of  a  given  population  from  the  sufficient 


statistics  for  the  normal  distribution  of  the  transformed  data  where  the 
logarithm  of  the  observed  values  is  taken  as  the  transformation  (the 
distribution  of  the  logs  is  known  to  be  normal). 

The  estimates  for  the  mean  and  variance  are  given  in  terms  of  an 
infinite  series  (actually  a  Bessel  function)  with  ;  ' ow  convergence  ex¬ 
cept  for  small  values  of  parameters  involved  ,und  these  are  modified  for 
better  arithmetic  computation. 

The  efficiency  of  estimations  without  considering  the  above  type 
transformation  is  also  discussed. 

36.  Fisherv  R.  A.,  and  Tippett,  L.,  "Limiting  Forms  of  the  Frequency 
Distribution  of  the  Smallest  and  Largest  Number  in  a  Sample," 

Proceedings  of  the  Cumb.  Phil.  Soc.,  Vol  24,  192?. 

Constants  for  obtaining  the  first  four  moments  of  the  distribution 

of  the  largest  member  of  n  sample  from  a  normal  population  for  samples 

up  to  1000  are  given,  possiole  limiting  forms  of  such  distributions  in 

general,  are  discussed.  Many  tabular  and  graphic  illustrations  are 

presented. 

37.  Furstenberg,  H.,  and  Kesten,  H. ,  "Products  of  Random  Matrices," 

Annals  of  Math.  Stat. ,  Vol.  31*  I960. 

R.  Bellman,  "Limit  Theorems  for  Non-Commutative  Operations,"  Duke 

Math.  Journal,  195^i  considers  this  problem:  Let  X  ,  X  ,  X  ,  form 

a  stationary  stochastic  process  with  values  in  the  set  of  k  x  k  matrices. 

He  investigated  the  asymptotic  behavior  of  ny^  =  Xn  Xn  ^  .  .  .  X^‘  and 

showed  that  if  the  Xn  are  independent,  then  the  lim  n  |log(nY''')i, 

n — ►  00 

exists,  (i,  j  refers  to  the  i  j«i  entry  of  the  mafcrixj  Bellman  conjectured 
that  n  '/'’(log(nY1)i ,  j  -  E  jlog  C^1)  i ,  3  |  )  is  asymptotics]. ly  normally 

distributed.  This  paper  adds  more  evidence  that  Bellman's  conjecture  is  correct. 
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38.  Geary,  R.  C.,  "The  Frequency  Distribution  of  the  Quotient  of  Two  Normal 
Variables,"  Journal  of  the  Royal  Stat.  Soc.,  Vol.  93 1  1930. 

Geary  supplies  an  approximation  of  q  -  x/y  which  is  widely  used.  He 

shows  that  for  X  and  Y  normally  distributed  with  respective  means  X  and  Y 

and  Q  defined  as  x./y#the  statistic 

Z  =  (q-QVCoC^2  q2  (■  O C.*  Q2)^ 

is  approximately  normally  distributed  where  o<^(i  =  1,  2)  represents 

|r  and  |r  . 
x  y 

39*  Geisser,  S.,  "The  Distribution  of  the  Ratio  of  Certain  Quadratic 
Forms  in  Time  Series,"  Annals  of  Math.  Stat.,  Vol.  28,  1957. 

In  testing  hypotheses  that  successive  members  of  a  series  of  ob¬ 
servations  are  serially  correlated, a  number  of  statistics  have  been 
proposed.  Durbin  and  Watson  gave  exact  distributions  for  some  of  these 
when  slightly  modified.  This  paper  extends  this  work  for  a  non-null  case 
of  two  of  these  modified  statistics  and  gives  a  simple  expression  for  the 
moments  of  another. 

40.  Goodspeed,  F.  M.,  "The  Relation  Between  Functions  Satisfying  A 
Certain  Integral  Equation  and  the  General  Watson  Transform" 

Canadian  Journal  of  ,Math. „  Vol.  2,  1950* 

Goodspeed  investigates  the  integral  equation 


r  00 

J  F(x)  F(u,x)dx  =  y-~ 


discussed  in  Chapter  VI.  In  particular,  a  class  of  functions  F(x) 
satisfying  the  equation  is  sought. 
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4l,  Gordon,  R.  D.,  "The  Estimation  of  a  Quotient  When  the  Denominator 
Is  Normally  Distributed,"  Annals  of  Math.  Stat.,  Vol.  12,  1941. 


This  investigation  is  devoted  to  obtaining  an  estimate  of  the  means 
of  components  of  ratio  values  obtained  in  a  time  series.  The  time 
series  is  of  the  form  x/y.  Of  corresponding  interest  in  the  time 
series  are  mean  values  x  and  y,  (true  value,  b/a) .  The  search  is  for 
an  estimate  of  (b/a)  which  has  the  property  that  it  may  be  numerically 
integrated.  Gordon  defines  a  function  Y(x)  so  that 

E  jy  >  y  (x)|  =  E(y)  •  E  {  7(x)  [  =  (b)  ^  =  b/a. 

The  function  derived  by  the  author  is  a  funotion  of  Mills’  ratio. 


42.  Grab,  E.  L. ,  "Tables  of  the  Expected  Value  of  1/X  for  Positive 
Bernoulli  and  Poisson  Variables,"  Journal  of  the  American  Stat. 
Assoc. .  Vol.  49,  195^ • 
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Tables  of  the  probability  function  of  (T  E  (;|  ,  where  X  is  a 
non -zero  random  variable  described  by  the  positive  Bernoulli  and  Poisson 
frequency  functions^re  presented.  Methods  of  solutions  are  presented 
which  are  more  efficient  than  the  factorial  series  method  first  pre¬ 
sented  by  Stephen [ll|  in  1945. 


43.  Greenwood,  J.  A.,  and  Greville,  T. ,  "On  the  Probability  of  Attaining 
a  Given  Standard  Deviation  Ratio  in  an  Infinite  Series  of  Trials," 
Annals  of  Math.  Stat.,  Vol.  10,  1958. 

Investigations  are  directed  toward  finding  the  probability  that  the 


ratio  of  the  deviation  from  the  expected  number  of  successes  in  a 

Binomial  experiment  of  n  trials,  (x^  -  np),to  the  standard  deviation 

y  n p q  (recomputed  after  each  trial)  will  exceed  some  positive  number  k. 

The  authors  prove  that  if  t^  =(xn~np)^^pq,  then  for  some  n  the  probability 

t  >  k  i8  unity, 
n  — 


44.  Greville,  T.  See  Greenwood  [43]  . 

45.  Gumbel,  E.  J.  See  Keeney  [64]  . 

46.  Gumbel,  E.  J,,  MThe  Geometric  Range  for  Distributions  of  Cauchy's 
Type,"  Annals  of  Math.  Stat. ,  Vol.  21,  1950« 

The  author  considers  the  distribution  of  the  reciprocal  of  the 
geometric  range  of  large  samples  drawn  from  population  characterized 
by  the  Cauchy  (or  nearly  Cauchy)  distribution.  Brief  tables  and 
graphs  of  the  probability  functions  are  presented. 

47.  Gumbel,  E.  J.,  and  Herbachi  L.  H. ,  The  Exact  Distribution  of  the 
Extremal  Quotient,"  Annals  of  Math.  Stat.,  Vol.  22,  1951. 

A  method  of.  obtaining  the  exact  distribution  of  the  quotient  ^  , 
of  the  extreme  values  found  in  a  sample  of  n  fc  2  observations  taken 
from  the  same  distribution  is  presented.  This  is  an  extension  of  the 
author's  first  published  work  (1950)  dealing  with  the  asymptotio  distri¬ 
bution  of  the  extremal  quotient.  The  authors,  in  this  paper,  consider 
the  Laplace,  Exponential,  Gamma,  Normal  and  Cauchy  distributions  as 
components  of  the  quotient  and  present  brief  graphs  of  the  prob¬ 
ability  function  associated  with  the  sample  size  n. 

48.  Gurland,  J.,  "Inversion  Formulae  for  the  Distribution  of  Ratios," 
Annals  of  Math.  Stat. ,  Vol.  19,  1948. 

Gurlaid  gives  theorems  and  proofs  that  the  distribution  function 

G(x)  of  (a,  x.  +  a_  x,  +  .  .  ,  +  a  x  )/(b.  x,  +  b„  x_  +  .  .  .  +  b  x  ) 
1122  n  n  1  1  2  2  nn 

is  obtainable  if  the  characteristic  function  of  the  joint  density 

function  (x^,  X£,  •  •  •  ,  xn)  is  known.  He  also  presents  inversion 

formulae  for  G(x)  and  shows  the  density  function  of  G(x).  Certain  other 

ratios  of  quadratic  forms  are  investigated. 


49.  Hack,  H.  R.  B.,  "An  Empirical  Investigation  into  the  Distribution 
of  the  F~Ratio  in  Samples  from  Two  Non-Normal  Populations," 
Biometrika,  Vol.  45,  1958. 

Experimental  data  obtained  from  an  agricultural  experiment  are 
used  to  estimate  the  distribution  of  the  F-ratio  selected  from  two 
highly  skewed  populations. 

50.  Hagis,  P.  See  Krishnaiah  [?oJ  . 

51.  Haldane,  J.  B.  S.,  "Moments  of  the  Distribution  of  Powers  and 
Products  of  Normal  Variates,"  Biometrika,  Vol.  32,  1942. 

Various  distributions  are  found,  the  derivations  of  which  aro 

greatly  detailed.  They  are:  (1)  the  distribution  of  a  oube  of  a 

normal  (0,  1)  variable,  (2)  the  general  case  of  the  distribution  of  any 

power  of  a  normal  variable,  (3)  the  distribution  of  the  product  of  n 

independent  normal  variables,  (4)  the  special  case  of  the  product  of 

3  normal  independent  random  variables,  (5)  the  product  of  two 

correlated  normal  variables,  and  (6)  the  Galton-Macalister  distribution, 

52.  Hart,  B.  I.,  and  VonNeumann,  J.,  "Tabulation  of  the  Probability  of 
the  Ratio  of  the  Mean  Square  Successive  Difference  to  the  Variance," 
Annals  of  Math.  Stat. ,  Vol.  13,  1942. 

2  2 

VonNeumann  determined  the  distribution  of  S  /  s  ,  the  ratio  of  the 
mean  square  successive  difference  to  the  variance  estimate,  for  odd 
values  of  sample  si&e  n.  In  this  paper,  the  probability  function 
developed  is  evaluated  for  other  values  (specific)  of  n.  The  evaluation 
methods  are  dependent  upon  the  Incomplete  Beta  functions. 

53.  Hartley,  H.  0.,  and  Ross,  A.,  "Unbiased  Ratio  Estimators,"  Nature, 
Vol.  1?4,  1954. 

Ratio  estimation  used  in  sample  surveys  to  estimate  the  population 
mean  Y  of  a  variate  y  with  the  help  of  the  known  population  mean  X  of 


some  correlated  variable  X  auffer  from  the  defect  that  they  are  biased 

estimators  and  by  amounts  for  which  ohly  approximate  formulae  are  offered. 

Confined  here  to  simple  random  samples  of  n  pairs  y^,  x^  from  a 

population  of  N  pairs,  various  ratio  estimators  of  Y  can  be  formed  such 

as  y  =  X  ^  ,  y1  =  Xr  where  y,  x,  r  are  arithmetic  means  of  samples  of 
x 

yi 

y. ,  x.  and  — ,  respectively.  To  get  exact  formulas  for  the  biases  one 

1  1  X. 


obtain©  E(y)  -  Y  -  Coy  ,  xj  f  E(y)  ~  Y  -  cov  ^  J  ,  x j  which  give 
bounds  on  the  biases:-  |bias  in  y  |  ;£  «T  y/x  0“  ,  )  bias  in  y]£  6~ y/x  -  CT  . 

A  j 

These  have  previously  been  attained  for  large  sample  approximation. 


54.  Hartley,  H.  0.  See  David,  H.  A.  [26]  . 


55.  Haviland,  E.  K. ,  "On  the  Distribution  Function  of  the  Reciprocal 
of  a  Function  and  of  a  Function  Reduced  Mod.  1,"  American  Journal 
of  Math.,  Vol.  65,  1941. 

Proofs  are  presented  which  show  that  under  suitable  conditions,  the 
functiunsl/x(  t )  and  x(t)  possess  asymptotic  distribution  functions  if 
x(t)  does.  x(t)  is  a  measur  real  function  defined  for  every  t. 

This  asymptotic  distribution  function  is  expressed  in  terms  of  the 
distribution  function  of  x(t)  for  certain  specified  cases. 

56.  Herbach,  L.  H.,  See  Gumbel,  [47]  . 

57.  Hess,  J.  (  See  Kish  [67]  . 

58.  Hirshfeld  ,  H.  0.,  "The  Distribution  of  the  Ratio  of  Covariance 
Estimates  in  Two  Samples  Drawn  from  Normal  Divariate  Populations," 
Biometrika,  Vol.  29,  1937. 

The  question  of  testing  whether  there  is  a  difference  between 
correlation  in  two  normal  populations  is  considered.  The  distribution 
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of  the  ratio  of  covariance  estimates  is  found  and  a  detailed  discussion 


i 

\ 
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of  the  appropriateness  of  its  use  in  this  test  is  presented.  The 
development  of  the  distribution  requires  frequent  use  of  the  incomplete 


Beta  function.  Some  agricultural  experiments  are  detailed. 


59*  Hogg,  R.  V.,  "On  Ratios  of  Certain  Algebraic  Forms,"  Annals  of  Math. 

Stat . ,  Vol.  22,  1951. 

Necessary  and  sufficient  conditions  for  the  stochastic  independence 
of  the  ratio  of  two  random  variables  and  its  denominator  are  given 
and  this  result  is  applied  to  special  linear  forms.  More  specifically 
let  x,  y  be  one  dimensional  random  variables  with  joint  density  function 
g(x,  y) .  Let  P(y  £  0)  =  0  and  assume  the  moment  generating  function 
M(  u,  t)  a  E  £exp  (ux+  t  y)  j  exists  for  -T^u,  t  <T,  T>0  ;  then  in 
order  that  y  and  r  =  -  be  stochastically  independent, it  is  necessary 
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and  sufficient  that 
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60.  Huntington,  E.  V. ,  "Frequency  Distributions  of  Products  and  Quotients," 
Armais  of  Math.  Stat.,  Vol.  10,  1939. 

Proofs  of  four  well  known  theorems  are  presented.  These  theorems  re¬ 
sult  in  a  mathematical  formulation  for  obtaining  the  distribution  of  a 
sum,  product,  quotient  and  difference  of  two  random  variables.  As  an 
example,  a  geometric  proof  is  given  to  show  that  the  distribution  of  the 
quotient  Z  =  x/y  is  given  by 


-oo  0 


fU,  y)  g(y)  y  dy 


where  x  and  y,(f(y)  t  0)  ,  are  independently  distributed  random  variables 

6l.  James,  Q.  S.,  "On  the  Accuracy  of  Weighted  Means  and  Ratios 
Biometrika,  Vol.  43,  1956. 


x^,  .  .  .  ,  are  k  quantities  derived  from  observations  which  are 

independent  and  normally  distributed  about  the  same  mean  y  but  with 
different  variances  <j^  ,  cr^  ,  *  *  *  * 

2 

are  known  but  estimates  s^  are  used.  A  confidence  limit  for  yi  is 

2 

required.  Let  w^  denote  1/  ,  the  reciprocal  -variance  6r  weight 

2 

of  ^and  let  z^  denote  l/(  ) .  Allow  w  a  and  z  =  Y.zi’  ** 

ia  reasonable,  in  the  absence  of  firm  knowledge  of  the  true  weights 

^  ,  _I1  ^ 
to  estimate  x  =  ^.w^^/w.  The  distribution  of  yu  =  w  (x  -  /i)  depends 

on  the  ratio  of  w^  .  The  desired  result  is  then  to  find  a  function 

/*  i  wp  \\ 

/t  I  ,  ”,  •  •  .  i  ^  J  such  that 

Pr  jjji  |  £  ji(w)j  =  P. 

This  function  is  found  and  tabulated  for  the  case  k  =  2. 

62.  Johnson,  N.  L. ,  "Some  Notes  on  the  Investigation  of  Heterogeometry 
in  Interactions,"  Trabajos  De  Estadistica,  Vol.  13,  1962. 

In  this  paper  certain  important  distribution  theory  of  use  in 
the  Analysis-of- Variance  model  for  analysis  of  interactions  ia  developed 


63.  Kac,  M.,  "On  a  Characterization  of  the  Normal  Distribution," 

American  Journal  of  Math. ,  Vol.  6l,  1939. 

A  aimplo  characterization  of  the  normal  distribution  is  presented. 

The  characterization  is  bused  on  an  invariant  principle  and  admits  a 
physical  interpretation. 

64.  Keeney,  R.  D. ,  and  Qumbei ,  E.  J.,  "The  Extremal  Quotient,"  Annals 
of  Math.  Stat. ,  Vol.  21,  1950. 

After  defining  the  extremal  quotient  as  the  ratio  of  the  largest  to 
the  absolute  value  of  the  smallest  observation  in  a  sample,  the 
authors  investigate  the  analytical  properties  of  symraetrioal,  continuous 
and  unlimited  distributions  through  a  ratio  defined  as  a  ratio  of  two 
non-negative  variates  with  identical  distributions.  Among  the  im¬ 
portant  results  is  the  fact  that  the  logarithm  of  the  extremalquotiont 
is  symmetrically  distributed.  Asymptotic  distributions  of  the  extremal 
quotients  from  the  Cauchy  and  Exponential  distributions  are  exhibited. 

65.  Keeten,  H.  See  Furstenberg  [37]  . 

66.  King,  E.  P.  See  Lukaos  [82]  . 

67.  Kish,  L. ,  and  Hess,  I.,  vi0n  Variances  of  Ratios  and  Their  Differences 
in  Multi-Stage  Sampling,"  Journal  of  the  American  Stat.  Assoc,, 

Vol.  54,  1959. 

The  aim  of  this  presentation  is  to  discuss  the  computation  for  the 
variances  for  the  estimators  r  =  x/y  and  (r  -  r’)  where  the  random 
variables  y  and  x  are  sample  totals  for  two  variables  obtained  from  w 
multi-stage  design.  The  variate  x  often  represents  the  sample  size  which 
represents  the  simplest  case  of  r.  This  and  several  other  cases  are  con¬ 
sidered.  Several  useful  computational  forms  are  presented  for 
var(r-r')  -  var  (r)  +  var  (r')  -  2  Cov  (r, r'). 
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68.  Kotlarski,  I.,  On  Random  Variables  Whose  Quotient  Follows  the 

Cauohy  Law,  Colloquium  Mathematloum,  Vol.  7,  I960. 

The  question  of  describing  a  claso  ^  of  random  variables  possessing 
the  property  thut  the  ratio  of  any  two  independent  random  variables  from 
the  class  having  the  same  frequency  function  follows  the  Cauohy  distribution 
is  considered.  A  method  is  given  for  constructing  an  arbitrary  number 
of  random  variables  belonging  to  along  with  the  necessary  and  sufficient 
conditions  for  belonging  to^f  .  Mellin  transforms  are  used  throughout 
the  analysis.  It  is  shown  that  the  random  variables  whose  frequency 
functions  are  determined  by  the  Mellin  transforms  which  are  solutions 
to  the  functional  equation  h(z)  h(-z)  =  - - — 


(-1*R  (z)*l) 

*  e  » 


COS  2 
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in  the  class  of  Mellin  transforms  of  symmetrical  random  variables  belong 
to  the  class  ^  . 

69.  Krishnaiah,  P.  R.,  P.  Hagis,  and  L.  Steinberg,  The  Bivariate  Chi 
Distribution,  Technical  Report  #3,  Applied  Math.  Dept.  Remington 

Rand  Univuc,  Philadelphia,  1961, 

70.  Krishnaiah,  P.  R.,  Hagis, P. ,  and  Steinberg,  L,  "a  Note  on  the 
Bivariate  Chi  Distribution,  "  SIAM  Review,  Vol.  5»  1963. 

Several  properties  of  the  bivariate  chi  p.d.f.  are  discussed  with 

some  mention  of  possible  applications.  Moments  of  the  .joint  p.d.f.  are 


presented  and  the  distributions  of  the  sum  of  U  and  V  and  the  ratio 
Z  =  U/V  are  shown. 


71.  Kruskal,  W.,  "The  Monotonicity  of  the  Ratio  of  Two  Non-Central  t 
Density  Functions,"  Annals  of  Math.  Stat.,  Vol.  25,  195^. 


The  ratio  of  two  different  non-oentral  t  density  funotiona  with 
the  same  number  of  degrees  of  freedom  is  shown  to  be  strictly  monotone 
with  sense  depending  on  the  relative  values  of  the  two  non-central 
constants.  The  author  also  cites  several  statistical  applications  in 
which  the  ratio  of  two  non-central  t  density  functions  arise. 


72.  Kuh,  E. ,  and  Mayer,  J. ,  "Correlation  and  Regression  Estimates  When 

the  Data  are  Ratios,"  Econometrics,  Vol.  23,  1955* 

The  authors  examine  the  influence  of  rat,io  transformations  on 
correlation  and  regression  estimates.  After  a  discussion  of  the 
"spurious"  ratio  correlation  problem,  neoessary  and  suffioient  conditions 
are  deduced  for  the  correlation  between  two  series  with  a  oommon 
denominator  to  equal  the  partial  correlation  between  numerator  series 
with  the  deflating  variable's  influence  held  constant.  It  is  shown  that 
conditions  must  be  fulfilled  to  obtain  the  beat  linear  unbiased  least 
squares  estimates  when  the  data  are  in  ratio  form. 


73*  Kullback,  S. ,  "The  Distribution  Laws  of  the  Difference  and  Quotient 
of  Random  Variables  Independently  Distributed  in  Pearson  Type  III 
Laws,"  Annals  of  Math.  Stat.,  Vol,  7,  1936. 


The  distribution  of  the  difference  u  =  x-y  where 


f(x) 


-x  P-1 
e  x 

~rTp y~ 


f(y) 


e-y  y0-:1 

~nqT“ 


for 


0;x  ;oo 
Of  y  ;f:oo 


is  first  found.  The  distribution  function  of  the  quotient  w  a  x/y  is 
investigated  through  the  utilization  of  the  transformation  w  a  eu  where 
u  =  log  x  -  log  y.  The  distribution  function  P(u)  is  found  and  with 
appropriate  substitution  the  solution  of  the  distribution  function  of  the 


quotient  is  found. 


74.  Laderman,  J.,  "The  Distribution  of  Student's  Ratio  for  Samples  of 
Two  Items  Drawn  from  Non-Normal  Universes,"  Annals  of  Math,  Stat., 
Vol.  10,  1939. 

The  formal  expression  for  the  distribution  of  student's  t  is 
derived  from  samples  of  two  drawn  from  any  population  having  a  con¬ 
tinuous  frequency  function.  A  geometric  method  similar  to  that  used 
by  Rider  is  used. 


75*  Laha,  R.  G.,  "An  Example  of  a  Non-Normal  Distribution  where  the 
Quotient  Follows  the  Cauchy  Law,"  Proceedings  of  the  National 
Academy  of  Science,  Vol.  44,  1958. 


For  x,  y,  independently  and  identically  distributed  random  variables, 


(the  density  function  of  x  being  f(x) 


-00  <.  x  <■  00  , ) , 


it  is  shown  that  the  ratio  x/y  follows  a  Cauchy  distribution  thus 
proving  untrue  a  previously  well  known  oonjeoture  that  this  was  a  unique 
property  of  normal  variates. 


76,  Laha,  R.  G. ,  "On  the  Laws  of  Cauchy  and  Gauss,"  Annals  of  Math. 

Stat. ,  Vol.  30,  1959. 

It  is  well  known  that  x/y  follows  the  Cauchy  distribution  when  x 
2 

and  y  are  normal  (0,  <r  ),  The  question  arises,  "is  it  possible  to 
obtain  a  characterization  of  the  normal  distribution  by  this  property 
of  the  quotient?",  Laha  constructs  a  class  of  functions  which  have  the 
property  that  their  ratios  are  Cauchy  distributed  but  the  components 
are  non-normal.  He  also  presents  and  proves  several  important  lemmas 
concerning  characterising  distributions. 

77.  Laha,  R.  G.,  "On  a  Class  of  Distribution  Functions  where  the  Quotient 
Follows  the  Cauchy  Law,"  Transactions  of  the  American  Math.Soc., 

Vol.  93,  1959- 

The  author  again  concerns  himself  with  the  distribution  of  the  ratio 


of  identically  distributed  random  variables  with  a  common  distribution 
function  F(x).  If  the  ratio  follows  the  Cauchy  law,  is  F(x)  necessarily 
normal?  Laha  constructs  a  different  class  of  functions  in  this  article. 

?8.  Loo'ye,  M. ,  ’’Fundamental  Limit  Theorems  on  Probability  Theory,”  Annals 
of  Hath.  Stat.,  Vol.  21,  1950. 

The  fundamental  limit  theorems  on  probability  theory  are  classified 
into  two  groups.  One  group  deals  with  the  problems  of  limit  laws  of 
sequences  of  sums  of  random  variables,  the  other  deals  with  the  problem 
of  limits  of  random  variables  in  the  sense  of  almost  sure  convergence. 

This  article  is  a  detailed  development  of  the  significant  results  in 
each  class  as  they  were  chrono logically  achieved. 

79*  Link,  H.  F.,  "The  Sampling  Distribution  of  the  Ratio  of  Two  Ranges 
From  Independent  Samples,"  Annals  of  Math.  Stat.:  Vol.  21*  1950. 

Investigation  of  the  joint  sampling  distribution  of  the  ratio  of 

ranges,  w^/w^*  where  w^  -  (xq  -  x^)/<r  ,  are  made  for  certain  specified 

distributions.  Tables  of  values  of  Prob  (w^w^  ^  R)  are  presented.  The 
underlying  distributions  considered  are  1)  Rectangular  Distribution, 

2)  e~X,  and  5)  Normal  Distribution. 

80.  Link,  R.  F. ,  "A  Correction  To:  The  Sampling  Distribution  of  the 

Ratio  of  Two  Ranges  From  Independent  Samples,"  Annals  of  Math.  Stat. , 
Vol.  23,  1952. 

A  correction  showing  that  the  probabilities  computed  in  the  reference 
above  are  correct  only  when  R  4  1  is  shown.  P.  R.  Rider  ("The  Distri¬ 
bution  of  the  Quotients  of  Ranges  in  Samples  from  a  Rectangular  Population," 
Journal.  American  Stat.  Assoc.,  Vol.  '45*  1951)  gives  the  correct  density 


function  for  R  >  1 
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81.  Lukaca,  E. ,  and  Szfisz,  0,,  "Certain  Fourier  Transforms  of  Distri¬ 
butions,"  Canadian  Journal  of  Mathematics,  Vol.  3»  1931 • 

The  question  of  when  a  given  function  can  be  the  characteristic 
function  of  a  probability  distribution  is  considered  for  a  restricted 
class  of  functions;  namely,  functions  which  are  reciprocals  of  poly¬ 
nomials,  with  the  view  .in  mind  of  doriving  conditions  which  are  easy  to 
apply.  The  following  necessary  conditions  are  derived; 

If  the  reciprocal  of  a  polynomial  without  multiple  roots  is  a 
characteristic  function,  then  the  following  two  conditions  are  satisfied: 

(1)  The  polynomial  has  no  real  roots.  Its  roots  are  either  all  on 
the  imaginary  axis  or  in  pairs  +  b  +  i  a  symmetric  with  regard  to 
this  axis. 

(2)  If  b  +  i  a  (a,  b  real;  a  /  0,  b  /  0)  io  a  root  of  the  poly¬ 
nomial  then  it  has  at  least  oile  root  1  o C  such  that  sgn  o£  =  sgn  a 
and  |  oL  I  <£.  |  a  I  . 

The  assumption  concerning  multiple  roots  is  used  in  deriving  (2^. 

82.  Lukacs,  E. ,  and  Kdng,  E.  P.,  "A  Property  of  the  Normal  Distribution," 
Annals  of  Math.  Stat.,  Vol.  25,  195^* 

The  authors  prove  the  theorem  under  Baau  [lO]  imposing  slightly  less 

restricted  conditions.  Their  conditions  require  only  that  the  random 

variables  x.. ,  x_,  .  .  .  ,  x  may  or  mav  not  be  identically  distributed 
J.  c  n 

and  possess  finite  moments  of  order  n. 

83.  Mann,  H.  B.,  and  Whitney,  D.  R-,  "On  a  Test  of  Whether  One  Variable 
is  Larger  than  the  Other,"  Annals  of  Math.  Stat,,  Vol.  18,  19^7. 

Let  x  and  y  be  two  random  vui’iableh  with  continuous  cumulative  distri¬ 
bution  functions  F  and  G.  A  statistic  u  depending  upon  the  relative  ranks 


of  the  x's  and  y's  is  proposed  in  this  paper  for  testing  the  hypothesis 
F  =  G.  The  conditions  under  which  f(u)  is  developed  are  stated  and  uses 
of  the  test  of  hypothesis  are  given. 

84.  Mayer,  J.  See  Kuh  [72]  . 


85.  Mauldon,  J.  G.,  "Characterizing  Properties  of  Statistical  Distri¬ 
butions.  "Journal  of  Math,  Oxford,  Vol.  7,  1956. 

Three  problems  dealing  with  the  general  theme  of  when  does  the 

distribution  of  a  sampling  statistic  determine  that  of  the  population 

are  discussed.  A  representative  problem  is  described  as  follows: 

Consider  X  and  Y  to  be  iridependently  distributed  variates  such 

that  r  =  X/Y  has  the  general  F  distribution. 
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dv,  O^v^oo;a,  b>0» 


Is  it  3till  true  that  X  and  Y  need  not  he  chi-square  variates?. 

86.  Merrill,  A.  S.,  "Frequency  Distribution  of  an  Index  when  Both  Com¬ 
ponents  !'■  How  the  Normal  Law,"  Biometrika,  Vol.  20,  1928. 

One  of  the  earliest  papers  to  investigate  the  deviation  of 

error  in  the  normal  distribution  approximation  lor  the  ratio  of  two 

random  variables.  In  this  paper,  the  author  investigates  the  ratio  of 

the  form  (y  +  y)/(x  +  x)  and,  through  graphic  means,  demonstrates  that 

when  correlation  is  high  between  x  and  y  and  the  coefficients  of 

variation  large,  there  is  a  considerable  deviation  from  normality. 

87.  Mickey,  M.  R.,  "Some  Finite  Population  Unbiased  Ratio  and  Re¬ 
gression  Estimators,"  Journal  of  the  American  Stat,  Assoc.,  Vol.  54, 

1959. 

A  class  of  ratio  and  regression  type  estimators  is  given  such  that 


the  estimators  are  unbiased  for  random  sampling,  without  replacement, 
from  a  finite  population.  Non-negative,  unbiased  estimators  of 
population  variance  are  provided  for  a  subclass  of  these.  Similar  re¬ 
sults  are  given  for  the  case  of  the  generalized  procedures  of  sampling 
without  replacement.  Efficiency  is  compared  with  comparable  sample 
selection  and  estimation  methods  for  this  case. 

88.  Mills,  J.  P. ,  "Table  of  the  Ratio:  Area  of  Bounding  Ordinate  for 
any  Portion  of  the  Normal  Curve,"  Bioaetrlka,  Vol,  18,  1926. 

The  tabular  values  i’or  .00  £x  10.00  (in  increments  of  .01)  of 


x 


and  a  discussion  of  its  derivation  and  many  possible  uses  are  presented. 

89.  Moore,  P.  G. ,  "Transformations  to  Normality  Using  Fractional  Powers 
of  the  Variable,"  Journal  of  the  American  Stat.  Assoc.,  Vol. 52,  1951* 

An  examination  of  the  properties  of  a  certain  class  of  trans¬ 
formations  is  made  under  the  assumption  that  they  are  designed  to 
transform  the  variate  X  into  some  form  of  the  normal  distribution.  The 
class  of  transformations  is  that  of 

x  ,  0  <  r  <  1  and  non-negative  X  . 

Some  consideration  is  given  to  a  Type  III  variate  X,  in  respect  to  the 
use  of  (X  -  Ot£)r  where  is  alao  a  random  va'riable  and  i3  distributed 
uniformly  on  the  interval  £0,  l]  . 

90.  Muity,  V.  N.,  "The  Distribution  of  the  Quotient  of  Maximum  Values 
in  Samples  From  a  Rectangular  Distribution,"  Journal  of  the 
American  Stat.  Assoc.,  Vol.  50,  1955- 

When  there  is  a  priori  knowledge  that  two  samples  have  been  drawn 
from  rectangular  populations  with  the  same  lower  hounds,  the  hypothesis 


that  both  samples  have  been  drawn  from  the  same  population  can  be 
tested  by  means  of  the  quotient  of  maximum  values.  The  distribution 
of  this  statistic  is  derived,  and  its  properties  studied.  Explicit 
expressions  for  the  power  function  of  the  test  are  given,  and  a 
table  of  %  values  of  the  quotient  is  presented  for  sample  sizes  up 
bo  ten. 

91.  Nair,  K.  R. ,  "The  Distribution  of  the  Extreme  Deviate  from  a 
Sample  Mean  and  Its  Standardized  Form,"  Biometrika,  Vol.35i  19^8. 

Denote  .....  x  as  a  random  sample  of  n  observations  drawn 

from  any  statistical  universe  so  that  the  observations  are  ranked  in 

ascending  order  of  magnitude.  The  author  first  summarizes  the  efforts 

of  other  studies  of  the  distribution  of  xy.  In  this  paper  attention 

is  concentrated  on  McKay's  method  of  solution  [Biometrika.  Vol.2?, 

p.  466]  .  It  is  shown  that  the  distribution  of  u  -  (x^-x^/  cr  ) , 

(termea  the  MoKay  statistic),  can  be  found  in  a  more  direct  way  than 

that  utilized  by  McKay.  The  distribution  of  u  is  reduced  to  certain 

integrals  termed  G-functiona.  Tabular  values  of  f(u)  are  presented. 

92,  Newman  ,  D. ,  "The  Distribution  of  the  Range  in  Samples  From  a 
Normal  Population  in  Terms  of  an  Independent  Estimate  of  the 
Standard  Deviation,"  Biometrika,  Vol.  31,  1939* 

The  author  makes  use  of  the  results  obtained  by  K.  Pearson  in 

computing  probability  levels  for  w/<T,  where  w  =  the  sample  range  and 

<T  =  population  variance ,to  determine  appropriate  corresponding  levels 

of  w/e  where  s  is  the  independent  estimate  of  CP  based  on  the  sample. 

The  distribution  of  q  =  w/s  is  obtained  and  factors  to  convert 

E(q)  to  E(w/cr  )  are  presented. 


93<  Nicholson  ,  C. ,  "A  Geometric  Analysis  of  the  Frequency  Distribution 
of  the  Ratio  Between  Two  Variables,”  Biometrika,  Vol.32,  1941. 


An  interpretation  of  the  geometry  of  the  distribution  of  (Y+Y)/(X+X) 

where  X,  Y  are  normal  (0,  <T ,  a~)  and  independent  is  presented.  The 

x  y 

case  of  dependency  is  given  some  consideration. 

94.  Noether,  G.  E. ,  "Two  Confidence  Intervals  for  the  Ratio  of  Two 
Probabilities  and  Some  Measures  of  Effectiveness,"  Journal  of  the 
American  Stat.  Assoo. ,  Vol.52,  1957* 

The  problem  of  comparing  a  suocess  rate  of  hypothetical  experi¬ 
mental  method  with  that  of  a  standard  method  is  attached  using  the 
statistic 


where  P^,  i  =  1,  2,  represent  the  success  ratio  of  the  old  and  new 
methods,  respectively.  This  is  a  departure  from  the  usually  method  of 
finding  a  confidence  limit  on  the  difference  A=  P2  “  Pi*  re” 

lative  advantages  and  disadvantages  are  discussed. 

93.  Olkin,  I.,  "Multivariate  Ratio  Estimation  for  Finite  Populations," 
Biometrika,  Vol.  45,  1958. 

This  paper  is  concerned  with  the  extension  of  ratio  estimation  to 
the  case  where  multi-auxiliary  variables  are  used  to  Increase  precision. 
The  following  model  is  presented: 


Population  Y^,  . 

*  *  ’  Yn  ; 

Y 

unknown 

Xll’  ' 

e 

*  ’  »  Xln5 

*1 

/  0  known,  R^  = 

e 

X  .  ,  . 

.  .  ,  X  *, 

• 

X 

/£  0  known.  R  = 

Pi 

pn 

P 

P 

and  the  (p+1)  x (p+1)  covariance  matrix  is  known. 
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The  proposed  estimator  of  Y  is  y  =  w,  r,  X.  +  .  „  .  +  w  r  X  where 

111  PPP 

w  =  (w. ,  .  .  .  ,  w  ),  X]w.  =1  is  a  weighting  function  and  r,  , 

p  i  X  Xi 

As  is  the  result  from  the  univariate  case,  y  is  shown  to  be 

biased  in  general  but  consistent.  The  Hartly-Ross  estimator  is 

generalized  to  yield  an  unbiased  estimator  of  Y. 

A  large  sample  approximation  is  given  for  the  mean,  variance  and 

_2 

mean  square  error  (to  0(n  ))  and  the  optimum  choice  of  the  weighting 

function  (ao  as  to  minimize  the  variance  of  y)  is  discussed. 

96.  O'Neill,  A.  See  Barton  [9]  • 

97.  Pearson,  K. ,  "On  the  Constants  of  Index  Distributions,"  Biometrika, 
Vol.  7,  1910. 

The  frequency  distribution  of  R  =  X/Y  is  investigated  by  calculating 
the  distribution  of  R  =  ZX  when  Z  =  l/Y.  The  known  properties  of  X  and 
Y  distributions  are  used  to  develop  expressions  for  the  momenta  of 
R  =  X/Y.  Pearson  applies  this  method  to  obtain  a  distribution  of  an 
Opsonic  index.  The  coefficients  of  variation,  means,  and  standard 
deviations  of  the  raw  data  are  known. 

■98.  Pearson,  E.  S.  See  David,  H.  A.  [26]  . 

99.  Pillai,  K.  C.  S.,  and  Ramachandran,  K.  V.,  "On  the  Distribution 
of  the  Ratio  of  the  i!t  Observation  in  an  Ordered  Sample  from  a 
Normal  Population  to  an  Independent  Estimate  of  the  Standard 
Deviation,"  Annals  of  Hath,  Stnt.,  Vol.  25,  1954. 

The  distribution  of  an  observation,  x^,  in  an  ordered  sample  of 
size  n  from  a  normal  population  with  zero  mean  and  unit  standard 
deviation  is  developed  as  a  series  of  Gamma  functions.  This  distri¬ 
bution,  in  turn,  is  utilized  to  find  the  distribution  of  q^  =  (x^/s) 


where  s  is  an  independent  estimate  of  the  standard  deviation  with  v 

degree  of  freedom.  In  a  similar  manner,  the  distribution  of  the 

"studentized"  maximum  modulus  u  -  lx  /si  is  obtained.  Uses  of  these 

n  1  n  1 

statistics  are  investigated. 

100.  Pothoff,  R.  F.  See  Roy  [ll3j  . 

101.  Das  Raj,  "On  Hill's  Ratio  for  Type  III  Populations,"  Annals  of 
Math.  Stat.,  Vol.  24,  19525. 

Several  authors  have  studied  the  Mills'  ratio,  i.e.,  the  ratio 
of  the  area  of  the  standardized  normal  curve  from  x  to  oo  and  the 
ordinate  at  x.  The  objective  of  this  note  is  to  establish  the  mono¬ 
tonic  character  of,  and  to  obtain  lower  and  upper  bounds  for,  the 
ratio  of  the  ordinate  of  the  standardized  Type  III  curve  at  x  and  the 
area  of  the  curve  from  x  to  oo  . 

102.  Romachandran,  K.  V  .  See  Pillai,  [99]  • 

103.  Reiter,  S.,  "Estimates  of  Bounded  Relative  Error  for  the  Ratio  of 
Variances  of  Normal  Distributions,"  Jouraal  of  the  American  Stat. 
Assoc.,  Vol.  51 •  1956. 

A  method  of  comparing  the  variability  of  two  competing  industrial 
processes  is  investigated.  In  this  paper,  the  process  whose  variability 

is  in  question  behaves  like  a  normally  distributed  random  variable  with 

2  2 
<J^  representing  the  variance  of  the  first  process  and  that  of  the 

2  2 

second.  The  ratio  cr^  /cn;  is  investigated  not  through  the  normal 
F  distribution,  but  through  methods  of  obtaining  a  bounded  relative 
error,  a  mini  max  estimate  as  suggested  by  Girshick  in  Theory  of  Games 
and  Statistical  Decision,  Wiley  and  Sons,  1954. 


104.  Rider,  P.  R. ,  "The  Distribution  of  the  Quotient  of  Ranges  in 
Samples  From  a  Rectangular  Distribution,"  Journal  of  the  American 

Stat.  Assoc.,  Vol.  22,  1951* 

The  distribution  of  the  quotient  of  the  ranges  of  two  independent 
random  samples  from  a  continuous  rectangular  distribution  is  investi¬ 
gated.  The  distribution  is  shown  to  be  independent  of  the  population 
range.  Uses  of  the  distribution  are  discussed.  The  special  case  of 
equal  sample  sizes  Is  investigated  and  four  brief  tables  of  probabilistic 
values  are  presented. 

105.  Rider,  P.  R. ,  "The  Distribution  of  the  Product  of  Ranges  in 
Samples  From  a  Rectangular  Population,"  Journal  of  the  American 
Stat.  Assoc. ,  Vol.  48,  1953* 

A  general  formula  is  given  for  the  product  of  the  ranges  of  two 
independent  samples  from  a  rectangular  population.  The  formula  does 
not  apply  to  cases  where  the  sample  sizes  are  the  same  or  if  they  differ 
by  unity.  Special  consideration  has  been  given  to  these  two  cases. 

106.  Rider,  P.  R. ,  "The  Distribution  of  the  Product  of  Maximum  Values 
in  Samples  from  a  Rectangular  Distribution,"  Journal  of  the 
American  Stat.  Assoc. ,  Vol.  50,  1955* 

This  note  is  to  supplement  the  article  by  Murty[90]  by  deriving 
the  distribution  of  the  product  of  maximum  values  in  random  samples 
from  a  rectangular  p.d.f.  The  case  of  the  product  of  K  maximum 
values  of  samples  of  equal  sizes  is  also  considered, 

107.  Rider,  P.  R.,  "Generalized  Cauchy  Distributions,"  Annals  of  Inst. 
Stat.  Math.,  Vol,  9,  195#  * 

This  paper  compares  the  asymptotic  variance  of  medians  in  samples 
of  size  n  from  any  distribution  (which  is  continuous  and  possesses  a  con¬ 
tinuous  first  derivative  in  the  vicinity  of  the  median)  with  exact 
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variances  from  distributions  of  the  following  type: 


f(x) 
where :  C 


It  is  seen  that  for  appropriate  values  of  k  and  h,  f(x)  becomes  the 
Cauchy  and  Student-Fisher  distributions. 


108.  Rider,  P.  R. ,  “Variance  of  the  Median  of  Samples  from  A  Cauchy 

Distribution, "  Journal  of  the  American  Stat.  Assoc..  Vol.55,  i960. 

The  exact  values  of  the  variances  of  the  medians  of  small  samples 

from  a  Cauchy  distribution  are  given.  The  tabular  values  have  been 

computed  from  the  integral  expression  representing  the  frequency 

function  of  the  median  of  a  sample  of ^2  K  +  lj  observations.  Values  of 

k  considered  are  integer  values  up  to  15. 


109.  Rider,  P.  H.,  "Expected  Values  and  Standard  Deviations  of  the 
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nomial  distribution  whose  zero  class  is  deleted.  A  recurrence  relation 


is  established  using  the  form 


x  k  f(x)  . 


x=l 


A  closed  form  of  E 

m 


(l) 


(the  expected  value  of  l/x) 


is  obtained  and 


results  are  computed  i--.8'.rjcally 
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u  =  (e  x  +  f)/(g  x  +  h) 


is  continuous  -1  £  x  <  1. 
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equation  rj- 


££---Jg 


(gx+h) 


2  * 


111.  Rietz,  H.  L.  "On  the  Frequency  Distribution  of  Certain  Ratios," 
Annals  of  Math.  Stat.,  Vol.  7*  1936. 

A  general  resume"'  of  the  published  articles  of  the  early  investi¬ 
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where  x,  n  are  integers,  x  ^  n, 

are  found  through  the  application  of  factorial  series.  Other  distri¬ 
butions  such  as  the  positive  Hypergeometrio  and  Poisson  are  considered. 
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xk 
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mean  square  successive  difference  S  "  be 
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.her.  .ia  -  >3i  =  |  *i><X;joC  V  (il  i  ‘  1 . n)’  ,h8n 
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a  sample  of  N  items  from  a  n-variate  normal  population  is  taken. 

N 

*i=  |  ^  is  the  sample  mean  of  the  iti  variate, the  value 

°c=  1 
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